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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Si)  units  as  follows: 


Multiply 


inches 

feet 

pounds  (mass)  per 
cubic  foot 

pounds  (force) 

pounds  (force)  per 
square  inch 

pounds  (force)  per 
square  foot 

feet  per  second 

degrees  (angle) 


JJL 


2.54 

0.3048 

16.01846 

4.448222 

6894.757 

47.88026 

0.3048 

0.01745379 


To  Obtain 


centimetres 

metres 

kilograms  per  cubic 
metre 

newtons 

pascals 

pascals 

metres  per  second 
radians 
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RESPONSE  AND  STABILITY  OF  EARTH  DAMS  DURING 


STRONG  EARTHQUAKES 

PART  I:  INTRODUCTION 

Background 

1.  With  the  failure  of  the  Lower  San  Fernando  Dam  during  the  San 
Fernando  earthquake  of  9 February  19719  the  attention  is  directed  towards 
the  problem  of  rational  design  of  earth  and  rockfill  dams  against  strong 
earthquake  ground  motion.  With  the  aid  of  a high-speed  computer  and  a 
very  sophisticated  method  of  finite  element  analysis  using  the  nonlinear 
properties  of  the  materials  at  the  time  of  failure  and  knowing  the 
ground  motion  characteristics  of  the  San  Fernando  earthquake,  it  is 
possible  to  back-analyze  and  give  the  cause  of  this  failure.  Seed  et 
al.  (1973).  But  for  the  designer,  who  has  yet  to  fix  a dam  section, 
to  satisfy  himself  about  the  soil  properties  obtained  from  laboratory 
tests  and  assume  the  ground  motion  characteristics  of  a future  event, 
the  use  of  the  finite  element  method  is  not  economical.  Moreover, 
the  results  obtained  by  a sophisticated  method  with  doubtful  input  data 
are  equally  doubtful.  Therefore,  the  use  of  a sophisticated  method 
with  doubtful  input  cannot  be  justified  unless,  of  course,  the  method 
is  equally  costly  and  can  be  used  with  equal  ease — not  the  case  with 
the  finite  element  method.  On  the  other  hand,  simplified  methods,  such 
as  those  based  on  linear  analysis  of  response  and  limit  equilibrium 
analysis  of  stability,  have  produced  results  comparable  to  those  of 
the  finite  element  method.  What  is  needed,  therefore,  is  further 
refinement  of  the  simplified  approach  without  adding  extra  cost  or 
time  to  the  designer  but  with  a more  extensive  use  of  information  about 
the  process  of  analysis.  Thus,  within  the  bounds  of  present-day  knowl- 
edge regarding  the  soil  parameters  and  the  ground  motion  input,  a 
simplified  method  of  approach  to  the  design  of  earth  da™s  is  desirable. 


h 


Purpose 


2.  The  purpose  is  to  provide  a comprehensive  method  for  such  an 
analysis.  In  order  to  analyze  the  stability  of  a dam  during  an  earth- 
quake, the  following  information  is  needed: 

ji.  The  inertia  forces  that  are  being  generated  in  the  dam 
during  an  earthquake. 

b.  The  resistance  of  the  dam  against  these  forces,  along 
with  the  preexisting  static  forces. 

£.  The  possible  consequences  when  the  resistance  is  not 
sufficient  to  withstand  these  forces  temporarily. 

These  three  factors  are  discussed  in  Parts  II,  III,  and  IV,  respectively. 

In  Part  V,  the  proposed  method  is  used  to  study  the  effect  of  the  San 

Fernando  earthquake  on  the  Lopez  Dam  in  California,  United  States  of 

America. 
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PART  II:  INERTIA  FORCES  AND  SEISMIC  COEFFICIENTS 
IN  EARTH  DAMS 

Inertia  Forces 

t 

3.  Earth  dams  are  not  absolutely  rigid,  and  when  excited  into 
oscillations  by  strong  earthquakes,  they  will  respond  in  a manner  that 
will  be  dictated  by  (a)  the  geometry  of  the  dam,  (b)  the  material 
properties  of  the  dam  and  its  foundation,  and  (c)  the  nature  of  the 
ground  movements.  Depending  on  circumstances , the  ground  movements  may 
cause  accelerations  to  develop  in  the  body  of  the  dam  with  associated 
fluctuating  stresses,  which  may  be  larger  or  smaller  than  those  of  the 
ground . 

1*.  During  a strong  earthquake,  which  may  last  from  a few  seconds 
to  a few  minutes,  the  ground  accelerations  will  fluctuate  with  time  in 
magnitude  and  direction.  A dam  to  which  this  excitation  is  applied 
will  respond  in  a manner  that  is  determined  by  its  resilience  and  its 
capacity  to  dissipate  energy.  The  lower  its  capacity  to  dissipate 
energy,  the  greater  will  be  its  response.  Moreover,  the  dam  will  seek 
out  the  periods  of  ground  movements  and  respond  strongly  to  those  with 
which  it  can  resonate.  Consequently,  an  earthquake  cannot  be  specified 
by  an  acceleration  alone,  nor  can  a dam  be  assumed  to  be  rigid. 

5.  To  make  the  problem  amenable  to  analysis,  it  is  necessary  to 
make  assumptions  to  obtain  such  solutions  as  various  authors  have  done. 

For  instance,  Ambraseys  and  Sarma  (1967)  give  a table  of  these  solutions, 
which  is  reproduced  here  for  easy  reference  (Table  l).  Ambraseys  and 
Sarma  (1967)  and  Sarma  (1968)  deal  with  the  problem  of  the  dam  on  the 
rigid  foundation  and  present  design  curves.  Since  it  is  not  always 
possible  to  assume  the  foundation  to  be  rigid,  it  is  thought  better  to 
produce  Mother  set  of  curves  which  will  take  the  depth  and  the  shear 
modulus  of  the  foundation  into  consideration. 

6.  Thus,  for  the  sake  of  analysis,  the  following  assumptions  are  j 

made  regarding  the  geometric  configuration  of  the  dam.  We  define  the 

geometry  by  Figure  1: 
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Solutions  for  Dam  Response  Problem 


Medvedev  and  Sinitsym  (1965) 
Clough  and  Chopra  (1965) 

p 

Unpublished  reports.  G(y)  = s(y)  p. 


9^  Rock 


Figure  1.  Geometry  of  dam- foundation  system 

a.  A triangular  untruncated  wedge  is  resting  on  a layer 
underlain  by  base  rock. 

b.  The  length  of  the  dam  is  much  greater  than  its  height. 

c_.  The  dam  is  symmetrical  about  a vertical  axis. 

d.  The  slopes  are  flat. 

e_.  It  vibrates  in  shear  only. 

7.  Under  these  assumptions,  a one-dimensional  shear  beam  analysis 
can  be  applied.  Slight  deviations  from  these  assumptions  do  not  affect 
the  result.  For  example,  Ambraseys  (i960)  has  shown  that  when  the  top 
width  of  the  dam  is  less  than  10  percent  of  the  base  width  of  the  dam, 
the  wedge  behaves  as  an  untruncated  wedge.  The  same  author  has  also 
shown  that  when  the  length  of  the  dam  is  larger  than  four  times  the 
height,  a one-dimensional  solution  applies.  Two-dimensional  solutions 
by  Clough  and  Chopra  (1965)  and  Hatanaka  (1955)  have  shown  that  for 
slopes  flatter  than  1:1.5  a shear  beam  solution  applies. 
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8.  Regarding  the  properties  of  the  material,  we  assume  that  the 
material  for  both  the  dam  and  its  foundations  are  linearly  elastic, 
homogeneous,  and  isotropic  but  that  the  properties  of  the  dam  and  the 
foundations  are  different.  Also,  the  material  has  an  energy  dissipating 
capacity  or  damping. 

9.  For  small  values  of  inertia  forces,  the  dam  behaves  elasti- 
cally, but  the  stress-  "train  behaviour  of  the  material  over  the  stress 
range  in  question  is  not  elastic.  However,  with  proper  use  of  damping 
values,  an  average  linear  relationship  may  be  used.  The  elastic  proper- 
ties of  the  dam  are  most  conveniently  introduced  in  terms  of  the  shear 
wave  velocity  of  the  material.  For  zoned  dams,  when  more  than  one  type 
of  soil  is  used,  a weighted  average  value  may  be  used.  Elastic  proper- 
ties of  soil  generally  vary  with  depth,  but  the  actual  variation  is  not 
exactly  known;  therefore,  its  use  cannot  be  justified  in  a simplified 
solution. 

Energy  dissipating  mechanisms 

10.  The  various  kinds  of  energy  dissipating  mechanisms  present  in 
the  dam  foundation  system  are  as  follows: 


a.  Viscous  damping.  In  this  case,  energy  loss  is  propor- 
tional to  the  velocity.  With  pure  viscous  damping  in 
elastic  material,  there  is  no  residual  displacement  at 
the  end  of  vibrations.  In  the  vibration  problems  for 
analysis,  this  is  the  most  convenient  form  of  damping. 

b.  Hysteretic  damping  due  to  nonlinearity.  This  is  the 
energy  loss  due  to  nonrecoverable  work  done  when  the 
material  is  strained  into  an  inelastic  zone.  In  this 
case,  there  may  be  permanent  deformations  of  the  struc- 
ture at  the  end  of  vibrations.  This  kind  of  damping 
varies  with  the  level  of  strain. 

c_.  Radiation  damping.  This  is  the  energy  loss  due  to  non- 
rigidity of  the  base  rock  when  some  of  the  vibrational 
energy  is  filtered  back  into  the  base  rock.  This  depends 
on  the  relative  properties  of  the  structure  and  the  base 
rock. 

In  reality,  all  three  energy-dissipating  mechanisms  are  present.  Since 
viscous  damping  provides  an  easy  method  of  solution,  we  lump  all  of 
these  into  one  equivalent  viscous  damping  term.  As  a result,  the  dis- 
placement that  is  computed  is  much  less  them  the  real  one.  On  the  other 
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hand,  the  stresses  and,  therefore,  the  inertia  forces  that  are  com- 
puted with  equivalent  viscous  damping  become  realistic.  Even  though 
the  pure  viscous  damping  is  a frequency-dependent  characteristic,  in 
vibration  problems  this  damping  is  assumed  to  be  constant  for  all 
frequencies. 

11.  There  is  no  direct  way  to  estimate  the  equivalent  viscous 
damping.  The  only  possible  way  is  vibrating  real  structures  and  re- 
calculating the  damping  from  the  response  of  the  structure.  Alterna- 
tively, it  may  be  possible  to  separate  the  radiation  damping  from  the 
other  two.  In  this  case,  the  radiation  is  considered  to  be  part  of  the 
problem,  but  the  solution  to  this  problem  does  not  exist  at  present. 

12.  The  nature  of  ground  movements  depends  on  various  factors, 
such  as  the  source  mechanism,  the  distance  of  the  site  from  the  source, 
and  the  geological  conditions  of  the  wave  path  from  the  source  to  the 
site.  The  radiation  of  energy  from  the  structure  to  the  base  rock  also 
affects  the  input  ground  motion. 

13.  The  earthquake  parameter  that  can  be  predicted  with  some 
confidence  is  the  maximum  ground  velocity.  Otherwise,  there  is  no  way 
of  predicting  a future  ground  motion.  One  possible  method  is  to  look 
through  the  library  of  existing  records  of  strong  motion  earthquakes 
and  find  the  one  that  was  produced  under  similar  conditions  of  source, 
distance,  and  geology.  An  alternative  is  to  select  a few  records  that 
have  different  frequency  characteristics  and  normalize  the  amplitudes 
to  obtain  the  expected  ground  velocity.  For  the  purpose  of  this 
analysis,  we  assume  that  the  earthquake  ground  motion  is  an  arbitrary 
disturbance. 

Derivation  of  governing 
equations  and  solutions 

Ik.  Let  us  consider  the  response  of  an  elastic  untruncated 
wedge,  resting  on  an  elastic  foundation  layer  underlain  by  a rigid 
base  (Figure  1).  We  consider  the  Vibrations  to  be  in  shear  in  one 
dimension  only,  i.e.,  perpendicular  to  the  y-axis  and  in  the  plane 
of  the  paper.  It  can  be  shown  that  the  equations  governing  the  vibra- 
tions of  the  dam  foundation  system  when  the  base  of  the  system  is 
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subjected  to  a ground  acceleration  g(t)*  are: 
a.  For  the  dam 


b.  For  the  foundation  layer 


where 

^Ug  = displacement  in  the  dam  and  the  layer,  respectively, 
measured  relative  to  the  base 

y = vertical  coordinate 

S^jSg  = shear  wave  velocity  in  the  dam  and  the  layer,  respectively 

••  •• 

U-^.Ug  = absolute  acceleration  of  a point  on  the  dam  and  the  layer, 
respectively 

c = damping  coefficient  for  both  dam  and  layer 

^l’^2  = velocity  °i"  a Point  in  the  dam  and  the  layer  relative  to 
the  rock,  respectively 

15.  These  assume  that  the  damping  coefficient  for  both  the  dam 
and  the  foundation  layer  is  the  same.  The  boundary  conditions  are 


u^h.t)  = u2(h,t) 

(3) 

Ug(H,t)  = 0 

(10 

3u 

57" = 0 at  y = 0 

(5) 

# For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  G). 
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M0(y)  • 


mJl(an) 

(<>%) 


sin 


cos 


(qan  IS) 


q = m(H-h)P2/(hP1) 


(12) 


(13) 


m = S1P1^2P2 


(1U) 


Sdn  = " IT  / ®(x)e  Xa)°n(t'T)  sin  [^U-t)]  dT 
n o 


w = u 
n on 


I 2” 

* / 1-X  = damped  and  undamped  circular  frequency, 

’ respectively,  of  the  dam  foundation  system 


“1  — 

u>o  = — an  = undamped  circular  frequency  of  the 

nth  mode  of  the  dam  foundation  system 


X = 2^ — = damping  as  a factor  of  critical, 
on  assumed  the  same  for  all  modes 


t = variable  of  integration 


(15) 


(16) 


(IT) 


(18) 


IT.  There  is  no  way  of  proving,  at  least  not  to  the  knowledge 
of  the  author,  that  the  solutions  u^(y,t)  and  u2(y,t)  obtained  by 
applying  the  Laplace  transform  technique  satisfy  the  governing  equa- 
tions. Therefore,  we  apply  the  limit  tests.  For  H=h  , the  solution 
u^(y,t)  converges  to  the  well-known  case  of  the  dam  on  the  rigid  base, 
which  is 


2J  (a  y/h) 

_ \ o nJ  s 

1 Z*  a J,  (a  ) dn 
, n 1 n 
n=l 


(19) 


where  a^  represents  n root  of  the  equation  J^a^)  = 0 . For  h=0  , 
the  solution  u^  converges  to  the  well-known  case  of  the  layer  on  a 
rigid  base,  which  is 
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“2  = 1 (-D-  • * S|te 

n=l 


(20) 


18.  There  is  another  indirect  test  that  can  be  applied  to  test 
the  solution,  which  is  based  upon  the  transient  part  of  the  solution. 
The  solution  obtained  through  the  Laplace  transform  technique  includes 
the  transient  part  of  the  solution.  For  the  undamped  solution,  the 
relative  displacements  u^  and  u^  should  equal  the  ground  displace- 
ment at  any  time  before  the  first  arrival  of  the  wave  to  the  point  in 
question.  Solution  u^  is  obtained  as  convolutions  of  two  functions 
F1  and  g , and  u2  is  obtained  as  convolutions  of  and  g , 

where  for  the  undamped  case 


v*  2J  (av/h) 

^(y.t)  = > r-2— 2 

L-l  anPo(q,m,n 


sin 


(Sl‘nt/h) 


Slan/h 


; h > y > 0 


(21) 


F2(y 


Z2M  (y) 

Z 

anP0(q.m,n 


sin 


(5l»nt/h) 


S a /h 
1 n 


; H >_  y >_  h 


(22) 


If  the  solution  is  true,  then 


F?(y,t)  = t for  0 <_  t <_  (H-y)/S„ 


(23) 


and 


F1(y,t)  = t for  0 < t < [(H-h)/S2  + (h-y)/sj 


(2U) 


19.  For  a set  of  typical  values  of  m , q , and  y , these  two 
functions  are  evaluated  numerically  for  t increasing  from  zero.  Func- 
tion F2  was  found  to  converge  rapidly  to  t during  the  period  con- 
cerned. For  values  of  y/h  > 0.2  , the  function  F^  also  converged  to 


lU 


t (Figure  2);  for  smaller  values  of  y/h  , particularly  for  y = 0 , 
the  convergence  is  very  poor,  and  it  is  assumed  that  it  will  converge 
to  t with  very  large  n , 

20.  If  the  two  solutions  are  correct,  then 


X 

I 


2J0(y/h) 

anPQ(q,m,n) 

2MQ(y) 

anPQ(q,m,n) 


!♦„  - 1 

1 


(25) 


(26) 


where  <t>  and  iji  are  the  n^*1  mode  shapes  for  the  dam  and  the  layer, 
n n 

respectively. 

21.  For  the  set  of  typical  values  of  m , q , and  y , the  con- 
vergence of  £<J>n  and  is  tested  and  found  to  be  not  uniform  at 

all  points  of  the  dam  (Figure  3).  However,  convergence  can  be  achieved 
with  very  large  numbers  of  modes  except  at  the  top  part  of  the  dam. 

22.  Double  differentiation  of  functions  u1  and  u2  now  yields 


U;L(abs)  = ^ + g = 1*^ 

(27) 

and 

u2(abs)  = u2  + g = 2>nSan 

(28) 

where  S 


2 ,2  2 
ui  - X o»  t 
n on 


an 


WN 


J g(T)e”*“on^'t_T^  sin  |un(t-T)]dT 


t .. 


+ 2Xu>  f g)t-t)  cos  fu>  (t-tjldT  (29) 

on  J L n J 


From  here  on,  we  shall  drop  the  term  (abs)  from  u-Jabs)  and  u (abs)  , 
••  •• 

and  simply  write  u^  and  u2  to  refer  to  absolute  accelerations,  but 


U1  ’ U2  and  “l  * and  ^2  Wil1  sti11  refer  to  relative  displacement 
and  velocity,  respectively. 

23.  In  view  of  the  nonuniform  convergence,  the  numerical  values 
of  u^  and  ug  are  obtained  in  the  following  way: 


and 


■i-  1 


<J>  S + 

Tn  an 


oc 


1 


. s 

n an 


(30) 


J 

■2 


ip  s 

n an 


I 

j+1 


ip  S 
n an 


(31) 


If  we  assume  that  the  response  acceleration  S (t)  of  a pendulum  of 

an 

period  T , such  that  T equals  0.05  sec  or  less,  is  equal  to  the 
ground  acceleration,  i.e.,  the  pendulum  behaves  as  a rigid  one,  then 


2 2 

\ - 2 *»s»  * 8u>  1 - 2 ♦» 

n=l  n+1 

U2  - 2 Van  * 8(t)  f1  - I *n) 
n=l  \ n=l  / 

and  j is  so  chosen  that 


(32) 


(33) 


Tj+^  0.05  sec 

2k.  Thus,  for  a given  ground  acceleration  g(t)  , the  complete 

response  time  history  can  be  computed  at  various  heights  of  the 

dam.  From  these  time  histories,  the  maximum  value  of  tt.  and  U„  

1 max  2 max 

can  be  determined  at  different  heights  of  the  dam  and  the  layer.  We  can 
determine  the  magnification  of  ground  acceleration  (gmax)  at  a point 

y = ah  given  by 
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1 max 


2 max 


Seismic  Coefficients 

25.  From  the  solutions  of  and  ti2  , it  is  seen  that  the 

accelerations  induced  in  the  dam  and  in  its  foundations  are  not  constant 
with  respect  to  height  or  in  time.  However,  most  of  the  stability 
analysis  techniques  prefer  to  have  the  accelerations  constant  with 
height.  According  to  Seed  and  Martin  (1966),  Ambraseys  and  Sarma  (1967), 
and  Sarma  (1968),  if  the  potential  sliding  surface  can  be  defined  in 
advance,  then  an  average  seismic  acceleration  can  be  computed  for  that 
surface  that  will  be  a function  of  the  ground  vibrations  and  periods  of 
the  dam. 

26.  The  "average  seismic  coefficient"  is  defined  as  that  factor 
which,  when  multiplied  by  the  weight  of  the  sliding  mass,  gives  the 
maximum  inertia  force  that  will  act  on  the  mass  during  a particular 
earthquake . Thus 


where 


K = average  seismic  coefficient 
Fmax  = maxiumum  inertia  force  on  the  potential  sliding  mass 
W = weight  of  the  sliding  mass 


For  convenience  of  comparison,  we  shall  express  the  seismic  coefficient 
in  terms  of  the  maximum  ground  acceleration  instead  of  the  acceleration 
due  to  gravity  as 

k.-Es-.Jjsb.  („) 


where  M represents  total  mass  of  the  sliding  body.  In  order  to 
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compute  K , therefore,  we  have  to  determine  the  total  inertia  force 
F and  the  total  mass  M of  the  sliding  body. 

27.  In  Figure  Ua,  consider  an  element  of  the  sliding  mass  of 
thickness  dy  and  width  b(y)  at  a depth  y from  the  crest.  The  mass 
of  the  elemental  slice  is 

dm  = pb(y)  dy  (36) 

where  p is  the  mass  density.  The  corresponding  absolute  acceleration 
at  the  depth  y at  an  instant  t is  U(y,t)  where  ti  = if  y h, 
and  u = if  y > h . Then  the  inertia  force  on  this  elemental  slice 
is 


dF  = U dm 

= pb(y)u(y,t)  dy  (37) 


Then,  the  total  mass  M of  the  sliding  body  is 

ah 

M = J pb(y)  dy  (38) 

o 

and  the  total  force  F is 

ah 

F = J pb(y)U(y,t)  dy  (39) 

o 

As  can  be  seen  from  Equation  39,  F is  a function  of  both  a and  t . 

We  can  evaluate  the  complete  time  history  of  the  function  F and  find 

the  absolute  maximum  value  of  it  at  some  instant  that  will  give  F , 

max 

for  a given  base  motion.  Thus,  knowing  F y , M , and  the  maximum 

value  of  the  ground  acceleration  g , we  can  determine  K , which 
0 max  a 

will  be  a function  of  a . For  convenience  of  computation,  we  may 
determine  a quantity 

A(a,t)  = | (UO) 
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and  then  compute 


A(a,t) 


K = 
a 


max 


(Ul) 


max 


One-parameter  sliding  wedge 

28.  In  order  to  compute  the  average  seismic  coefficient,  it  is 
essential  to  define  the  potential  sliding  surface  in  advance.  The  most 
convenient  is  the  one  suggested  by  Seed  and  Martin  (1966) , Ambraseys 
and  Sarma  (1967),  as  shown  in  Figure  Vb.  These  sliding  wedges  have  a 
horizontal  base,  and  the  slip  surfaces  pass  through  the  apex  of  the  dam. 
The  seismic  coefficients  computed  for  these  surfaces  become  independent 
of  the  base  width  of  the  sliding  wedge,  and  these  are  also  independent 
of  the  slopes  of  the  dam.  The  only  parameter  that  is  required  to  de- 
scribe such  a surface  is  a the  depth  of  the  base  of  the  sliding 
surface  from  the  crest.  As  can  be  seen,  these  surfaces  can  be  defined 
only  in  the  dam  and  the  foundation  layer  cannot  be  included. 

For  such  a surface. 


b(y)  = By 


where  B is  a constant  and 


Then 


A(a,t)  = 


«<y,t)  - Vir.t)  = £ ^ * g(t)  U - 

ah  j / J \ 

/ y L *nSan  + 8(t)  V1  * I Ody 

o nfl \ n=l  / 

ah 

/ ydy 


(U2) 


(1*3) 


J ah  ah  j ah 

t S&n  J y V ^ + B(t)  f Y ~ 8(t)  £ f y 4>n  dy 

n=l  o o n=l  o 


ah 

/ ydy 
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Then 


A(ct,t)  = V *1  (o)Sari  + g(t) 
4*.  n an 
n=l 


t1 ' l H 


(1*1*) 


where 


ah 


4>1  (a)  = — 
n'  <*h 


/ y+n  dy 


/ y dy 


(1*5) 


As  shown  in  the  Appendix  B 


U J ( aa  ) 

♦V»>  - -V— -■ 


aanPo(q,m,n) 


(i *6) 


Four-parameter  sliding  wedge 

29.  Because  of  the  inability  of  the  one-parameter  sliding  wedge 
to  include  a part  of  the  foundation  layer,  Ambraseys  and  Sarma  (1967) 
devised  the  four-parameter  sliding  wedge  as  shown  in  Figure  1+c.  In 
these  wedges,  the  four  parameters  are  the  depth  of  the  point  C from  the 
crest  (a),  the  slope  of  the  dam,  and  the  two  angles  0 and  \Ji.  It  was 
assumed  that  the  position  of  the  point  C is  fixed.  But  in  reality,  the 
position  of  the  point  C can  also  be  changed,  which  of  course  increases 
one  parameter  and  seems  unnecessary. 


30.  For  such  a surface,  assuming  = p 


ah 


F = J pb(y)U1  dy 


0 < a < 1 


(1*7) 


or 


and 


h ah 

F * f pb(y)U1  dy  + C pb(y)U2  dy  if  a 
M * f pb(y)  dy 


> 1 


(1*8) 


(1*9) 
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This  also  yields. 


A = T 4>U  (o)S  + g(t) 
a r n an 

n=l 


t - i 


where 


$1*  (a)  = 
n 


h qh 

J*b(y)^n  dy  + J bfy)^  dy 


for  a > 1 


ah 

J"  b(y)  dy 


See  Appendix  B for  «J>A  (a)  . 


Discussions  and  Design  Curves 


31.  The  single  parameter  by  which  a dam-layer  system  can  be 
represented  is  its  fundamental  period.  The  fundamental  period  and  the 
higher  modes  are  functions  of  four  factors  - S1  , h , m , q . These 
represent 


(a)  The  shear  wave  velocity  of  the  material  of  the  dam  = S.^ 

(b)  The  height  of  the  dam  = h . 

(c)  The  shear  wave  velocity  of  the  foundation 

pl^l 

layer  - m . 

P2 

(d)  The  height  of  the  layer  on  which  the  dam  is  situated 
h . 


= il!l 


m p. 


The  fundamental  period  of  the  dam-layer  system  is  then  given  by 


T = 2rr 

a!S! 

where  a^  is  the  1st  root  of  the  equation 

J (a  ) 

m tan(qa  ) = — — — 

n w 

32.  The  higher  roots  of  the  equation  give  the  higher  mode 
periods  of  the  system.  These  periods  are  given  by 


(56) 


(57) 


Tn  - T1*A 


(58) 


Figure  5 gives  the  values  of  a1  for  a combination  of  values  of  m and 
q . These  curves  can  be  used  to  determine  the  fundamental  period  of 
the  dam-layer  system.  Some  higher  roots  of  the  equation  are  given  in  a 
tabular  form  in  the  Appendix  C. 

33.  The  following  earthquake  acceleration  records  are  chosen  in 
order  to  determine  the  maximum  accelerations  at  various  heights  of  the 
dam.  These  records  were  chosen  mainly  because  of  their  differences  in 
frequency  characteristics.  The  ground  motions,  acceleration,  velocity. 
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Figure  5.  Fundamental  period  of  dam-foundation  system 


and  displacement  spectra  of  these  records  are  shown  in  Figures  D1  through 
D9  (Appendix  D).  An  example  (the  Port  Hueneme  record)  is  shown  in 
Figure  6.  The  records  axe 

a.  San  Fernando  Earthquake  1971  - Pacoima  Dam  3 components 

b.  Imperial  Valley  Earthquake  I9U0  - 

El-Centro  Record  2 components 

£.  Koyna  (India)  Earthquake  196?  - 

Koyna  Dam  Record  2 components 


Another 

later. 


d..  Parkfield  Earthquake  1966  - Cholame, 

Shandon  Array  No.  2 1 component 

£.  Port  Hueneme  Record  1957  1 component 

record,  which  is  used  for  the  study  of  the  Lopez  Dam  is  added 

This  record  is 

f\  San  Fernando  Earthquake  1971  - Deconvoluted  and 
modified  version. 
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Figure  6.  Ground  motion  records  and  spectra:  Port  Hueneme  record  of  1957  (sheet  1 of  2) 
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Figure  6 (sheet  2 of  2) 


34.  In  Figures  7 and  8 the  effects  of  m and  q are  shown.  It 
appears  that  for  m less  than  0.2,  there  may  be  some  convergence 
problems,  particularly  in  the  upper  part  of  the  dam,  when  the  predomi- 
nant period  of  the  earthquake  ground  motion  nearly  coincides  with  the 
fundamental  period  of  the  dam-layer  system.  For  values  of  m less  than 
0.7,  there  is  some  magnification  in  the  response  of  the  dam,  compared 
with  the  response  of  the  dam  on  the  rigid  base.  The  magnification  is 
largest  for  values  of  m about  0.25.  The  magnification  is  also  af- 
fected by  the  value  of  q/m  as  shown  in  Figure  8.  The  magnification 
is  different  at  different  heights  of  the  dam  and  at  different  periods 
of  the  dam-layer  system. 

35-  Figure  9 is  an  explanation  of  Figures  E-2  through  E-127  (Ap- 
pendix E)  that  show,  in  spectra  form,  the  maximum  response  of  dams  of 
various  fundamental  periods,  and  m and  q factors  subjected  to  the 
first  nine  of  the  earthquake  records  listed  above.  The  responses  are 
calculated  for  a damping  value  of  20  percent  critical  for  all  modes.  A 
maximum  of  20  modes  was  used.  The  values  of  m and  q that  are  chosen 
are 


m = 0 

0.5 

i.o 

q - 0 

0.125 

0.25 

0.375 

0.50 

0.75 

1.0 

0.25 

0.50 

0.75 

1.00 
1.50 

2.00 

36.  These  sets  of  q values  correspond  to  depth  of  layer  to 
height  of  dam  ratios  of 

= 1/4,  1/2,  3/4,  1,  1-1/2  and  2 for  P;L  * p2 

At  the  same  time,  the  seismic  coefficients  for  one-parameter  and  four- 
parameter  sliding  wedges  are  computed  and  shown  in  Figure  9 in  spectra 
form.  In  this  figure.  Figure  9a  gives  the  values  of  the  amplication 
factor  Ma  (the  maximum  acceleration  at  a point  divided  by  the  maximum 
ground  acceleration),  for  values  of  o = 0,  0.2,  0.4,  0.6,  0.8,  1.0,  1.2, 
1.4,  1.6,  1.8  if  the  depth  of  the  layer  is  greater  than,  or  equal  to 
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AMPLIFICATION  FACTOR  Ma  AMPLIFICATION  FACTOR,  M 


j 
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the  height  of  the  dam.  For  smaller  depths  of  layer,  the  number  of 
curves  is  curtailed  at  the  corresponding  point.  For  example,  q/m  =0.5 
represents  a layer  depth  which  is  half  the  height  of  the  dam.  Therefore, 
a = l.U  is  the  last  curve  as  the  total  depth  from  the  crest  of  the  dam 
to  the  base  rock  if  y = 1.5h  . 

37.  Figure  9b  gives  the  seismic  coefficients  for  one- 

parameter  sliding  wedges  (seis.  coef.  1)  for  values  of  a = 0.2,  0.4. 

0.6,  0.8,  and  1.0  . 

38.  Figure  9c  gives  the  seismic  coefficients  for  four- 

parameter  sliding  wedges  (seis.  coef.  4)  for  values  of  a = 0.2,  0.4, 

0.6,  0.8,  1.0,  1.2,  1.4,  1.6,  1.8,  2.0  when  the  layer  is  greater  than, 
or  equal  to,  the  height  of  the  dam.  For  smaller  layer  depths,  the 
number  of  curves  is  curtailed  at  the  corresponding  a level,  the  last 
curve  being  an  odd  one  representing  a wedge  touching  the  base  rock. 

For  example,  for  q/m  = 0.5  » the  next  to  last  curve  is  for  a = 1.4  , 
and  the  last  curve  is  for  a = 1.5  . The  other  parameters  for  the 
four- parameter  wedges  are  0 = 45  deg,  ip  = 80  deg,  slope  = 1:2.5  . 

39.  Figure  10  is  an  explanation  of  Figures  F-2  through  F-15  (Ap- 
pendix F)  that  show  the  envelopes  of  the  nine  earthquake  acceleration 
records.  These  envelopes  can  be  used  to  determine  the  maximum  values 
of  seismic  coefficient  when  the  frequency  characteristics  of  probable 
ground  motion  are  unknown. 


MAXIMUM  OF  9 RECORDS 


PART  III:  STABILITY  ANALYSIS 


Critical  Acceleration  of  Slip  Surfaces 

1+0.  In  Part  II,  the  inertia  forces  and  the  seismic  coefficients 
are  determined  for  the  dam  and  its  foundations.  The  problem  then  is  to 
determine  whether  the  dam  section  will  be  able  to  withstand  these  forces. 

1+1.  A very  convenient  term  to  describe  whether  a given  slip 
surface  will  be  able  to  withstand  these  forces  is  the  critical  accelera- 
tion. The  critical  acceleration  is  defined  as  that  horizontal  acceler- 
ation which  is  required  to  bring  a slip  surface  to  a state  of  limiting 
equilibrium,  and  this  term  is  therefore  analogous  to  the  factor  of 
safety.  If  the  critical  acceleration  for  the  slip  surface  is  larger 
than  the  corresponding  seismic  coefficient,  then  the  surface  has  a 
factor  of  safety  greater  than  one  and  is  therefore  safe.  In  contrast  to 
the  static  problem  of  determining  the  minimum  factor  of  safety  for  the 
entire  section,  a minimum  critical  acceleration  for  the  section  for  the 
earthquake  problem  is  not  meaningful.  As  shown  in  Part  II,  the  seismic 
coefficients  for  slip  surfaces  at  different  heights  are  different.  The 
surface  that  produces  the  minimum  critical  acceleration  could  be  safe 
against  the  corresponding  seismic  coefficient,  but  another  surface  with  a 
higher  critical  acceleration  may  not  be  safe,  particularly  in  the  upper 
part  of  the  dam.  Therefore,  it  will  be  necessary  to  determine  the 
critical  acceleration  for  slip  surfaces  at  various  heights  of  the  dam. 

1+2.  According  to  Sarma  and  Bhave  (197^)»  without  the  change  of 
strength  characteristics  and  change  of  pore  pressures,  the  critical 
acceleration  for  any  given  slip  surface  bears  a relationship  with  the 
factor  of  safety;  therefore,  the  critical  acceleration  can  be  used  as  a 
measure  of  the  static  factor  of  safety.  However,  any  change  in  these 
conditions,  such  as  the  increase  of  pore  pressures  due  to  earthquake 
loading,  removes  these  relationships.  It  has  been  shown  (Sarma,  1973) 
that  for  a given  dam  section  and  for  static  condition  or  pseudo-static 
condition  of  earthquake  loading  with  horizontal  accelerations,  the  de- 
termination of  the  critical  acceleration  is  easier  than  the  determination 


of  the  factor  of  safety.  The  advantages  of  using  a critical  accelera- 
tion solution  instead  of  a direct  factor  of  safety  solution  're  men- 
tioned in  Sarma  and  Bhave  (197*0. 

1*3.  In  this  analysis,  it  is  assumed  that  the  earthquake  accelera- 
tion is  constant  along  the  height.  It  is  known  that  the  earthquake 
acceleration  is  not  constant  along  the  height,  and  that  is  why  the 
seismic  coefficient  was  defined  in  Part  II.  Therefore,  an  analysis 
is  made  to  study  the  effect  of  the  acceleration  distribution  on  the 
critical  earthquake  load  (Valenzuela,  1975).  From  this  limited  study, 
it  is  found  that  the  critical  earthquake  load  for  a given  slip  surface 
is  minimum  when  the  acceleration  distribution  is  constant  and  is  not 
minimum  when  the  accelerations  are  increasing  towards  the  crest.  Thus, 
the  use  of  a constant  acceleration  with  the  seismic  coefficients  for  that 
surface  produces  conservative  results.  But  the  use  of  the  critical 
acceleration  with  constant  distribution  cannot  produce  the  critical 
surface,  whereas  the  use  of  varying  acceleration  can.  However,  the  use 
of  varying  acceleration  for  design  purposes  is  not  suggested,  for  the 
variation  is  not  known  in  advance. 

Existing  Methods  of  Stability  Analysis 

1+U.  It  is  accepted  that  the  pore  water  pressures  in  a dam  increase 
with  cyclic  loading  induced  by  earthquakes.  Since  increase  in  pore 
pressures  reduces  the  total  strength,  this  plays  an  important  part  in 
the  stability  of  the  dam  section.  But  before  going  into  the  solution 
of  the  problem  for  earthquake  loading,  it  is  necessary  to  examine  the 
problem  of  the  stability  analysis  both  for  static  and  earthquake  condi- 
tions. What  follows  below  is  an  assessment  of  this  problem.  From  this 
we  will  see  that  different  techniques  can  be  developed  that  will  be 
suitable  for  specific  types  of  problems.  For  example,  to  analyze  an 
existing  slip  surface  in  a dam  or  a slope,  the  technique  of  solution 
should  be  different  from  that  required  for  designing  a dam.  Since  in 
the  design  of  dams,  one  can  reject  a slip  surface  on  the  basis  that 
the  surface  is  not  kinematically  admissible,  the  existing  slip  surface 


must  produce  a kinematically  admissible  solution.  Thus,  the  assumptions 
may  be  varied  to  suit  the  problem.  But  whatever  the  technique,  the  so- 
lution is  based  on  the  limit  equilibrium  principle. 

1*5.  There  are  several  methods  of  stability  analysis  in  existence. 
All  of  these  have  certain  things  in  common: 

a.  These  methods  are  based  on  the  limit  equilibrium  principle 
which  defines  the  factor  of  safety  or  the  critical 
acceleration. 

b.  They  employ  the  Mohr-Coulomb  failure  criterion  with 
rigid-plastic  models. 

c.  All  of  these  apply  a method  of  slices. 

d.  They  are  based  on  assumed  slip  surfaces. 

These  methods  differ  in 

e.  The  shape  of  the  assumed  slip  surfaces. 

f.  The  handling  of  the  indeterminancy  of  the  problem  as 
shown  later. 

H6.  Before  proceeding  further,  we  need  to  define  a few  terms. 

a.  Simplified  solution.  A solution  which  does  not  satisfy 
all  the  static  equilibrium  conditions.  Assumptions  are 
made  to  obtain  the  solution  in  a simple  form.  Most  of 
the  available  solutions,  including  Bishop's  simplified 
method  (1955),  Kenney's  method  (1956),  and  even  Janbu's 
generalized  method  of  slices  (1957),  fall  into  this 
category. 

b.  Rigorous  solution.  A solution  which  satisfies  the 
complete  equilibrium  condition  of  statics.  The  implied 
forces  obtained  from  the  solution  may  not  be  acceptable. 

c.  Acceptable  solution.  A rigorous  solution  which  has  the 
added  criterion  of  the  implied  forces  satisfying  the 
criterion  of  acceptability.  The  basic  criteria  are: 
the  forces  must  not  violate  the  Mohr-Coulomb  failure 
law,  no  tension  is  implied,  and  the  direction  of  forces 
must  be  kinematically  admissible. 

Our  aim  is  to  obtain  an  acceptable  solution  with  the  least  possible 
effort. 

U7.  In  the  method  of  slices,  the  possible  failure  mass  is  divided 
into  n slices  (Figure  ll).  For  each  slice,  the  known  and  unknown  forces 
and  their  points  of  applications  are  shown.  For  n slices,  we  have  the 
following  (6n-2)  unknowns: 


n numbers  of  the  normal  force  N 
n numbers  of  the  shear  force  T 
(n-l)  numbers  of  the  body  force  E 

(n-l)  numbers  of  the  body  force  X 

(N-l)  numbers  of  the  point  of  application  of  E 

n numbers  of  the  point  of  application  of  N 

1 in  the  form  of  factor  of  safety  or  critical 
accelerations 


c — 
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Figure  11.  Stability  analysis  by  method  of  slices 

48.  From  the  static  equilibrium  conditions  of  each  slice,  we 
have  the  following  equations: 


£ Moment  = 0 
2 Vertical  forces  = 0 
£ Horizontal  forces  = 0 


The  Mohr-Coulomb  failure  criterion  gives  T = f(N)  ; hence,  for  n 
slices  we  have  4n  equations.  In  order  to  obtain  a rigorous  solution, 
therefore,  one  has  to  make  2n-2  assumptions.  If  more  assumptions  are 
made,  more  unknowns  must  be  introduced.  Now  let  us  look  into  the  various 
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available  solutions  and  see  how  these  assumptions  are  made. 

a.  Bishop's  simplified  method:  In  this  method,  the  assump- 
tions are  n numbers  of  points  of  application  of  N , 
and  (n-l)  numbers  of  the  magnitude  of  X (X  = 0).  Since 
one  more  assumption  is  made  than  required  and  no  unknown 
is  introduced,  one  known  condition  must  be  sacrificed. 

In  this  case,  the  condition  that  is  not  satisfied  is 
En+1  # 0 . For  circular  slip  surfaces  this  error  is 
small,  and  the  result  is  therefore  mostly  acceptable. 

b.  Janbu's  generalized  method  of  slices:  The  assumptions 

are  n numbers  of  points  of  application  of  N , and 

(n-l)  numbers  of  points  of  application  of  E . Again 

there  is  one  more  assumption  than  is  required.  In  this 

case,  the  condition  that  is  not  satisfied  is  M i 0 for 

n 

the  last  slice. 

c.  Kenney's  method:  The  same  assumption  as  Janbu's  is  made, 

and  the  equation  that  is  not  satisfied  is  ^ 0 • 

d.  Morgenstern  and  Price  method  (1966):  The  assumptions 
are  n numbers  of  point  of  application  of  N , and 
(n-l)  numbers  of  the  relationship  between  X and  E . 
Since  they  have  made  one  assumption  more  than  is  required, 
an  extra  unknown  X is  introduced.  This  method,  there- 
fore, satisfies  static  equilibrium  conditions  rigorously. 

e.  Sarma's  method  (1973):  The  assumptions  are  n numbers 
of  points  of  application  of  N , and  (n-l)  numbers  of 
relative  magnitudes  of  X . As  in  Morgenstern  and 
Price's  method,  an  extra  unknown  X is  introduced; 
therefore,  it  is  a rigorous  solution. 

1+9.  The  necessary  and  sufficient  number  of  independent  assumptions 
that  are  required  are  (2n-2).  If  one  can  make  these  assumptions  satis- 
factorily, a rigorous  solution  can  be  obtained.  Some  possible  combina- 
tions of  assumptions  that  one  can  make  are  given  below. 

a.  (n-l)  numbers  of  relationship  between  X and  E , and 
(n-l)  numbers  of  point  of  application  of  E or  N . A 
method  of  solution  baaed  on  these  assumptions  is  provided 
herein. 

b.  (n-l)  numbers  of  point  of  application  of  E , and  (n-l) 
numbers  of  point  of  application  of  N . This  means  that 
Janbu's  method  can  be  a rigorous  one  by  assuming  that  the 
point  of  application  of  the  last  N is  unknown  and  then 
determining  it  as  part  of  the  solution. 

c.  (n-l)  numbers  of  magnitudes  of  X , and  (n-l)  numbers  of 
point  of  application  of  E or  N . This  method  is  not 
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recommended,  as  the  solution  depends  on  the  absolute 
magnitude  of  X , which  may  be  very  difficult  to  guess. 

d.  n numbers  of  relative  magnitude  of  N and  (n-l)  numbers 
of  point  of  applications  of  E or  N introduce  an 
extra  unknown  in  the  form  of  X . Since  this  is  a very 
simple  solution,  it  is  also  provided  here. 

New  Methods  of  Stability  Analysis 

50.  It  may  be  possible  to  find  other  sets  of  assumptions,  but 
whether  any  one  of  these  assumptions  will  provide  acceptable  solutions 
is  debatable.  What  follows  below  are  two  new  methods  of  solutions  of 
stability  of  slopes. 

Method  1 

51.  Figure  12  shows  a slip  surface  in  a given  slope.  The  slip 
surface  is  formed  of  segments  of  n straight  lines.  We  subdivide  the 
sliding  mass  into  n slices,  but  the  slices  need  not  be  vertical  and 
the  edges  of  the  slices  need  not  be  parallel.  The  force  on  the  slice  i 
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is  shown  in  the  inset.  We  assume  that  under  the  action  of  the  horizontal! 
acceleration  kcg  , the  slip  surface  is  in  limiting  equilibrium. 

52.  From  the  vertical  and  horizontal  equilibrium  of  the  slice, 

N.  cos  a.  + T.  sin  a.  = W.  + cos  - X.  cos  6. 

1 l i i l l+l  l+l  l l 

- Ei+1  sin  6i+1  + E^  sin  (59) 

where 

a = the  slope  of  the  base  of  the  slice  with  respect  to  the  hori- 
zontal 

6 = slope  of  the  side  of  the  slice  with  respect  to  the  vertical 

T.  cos  a.  - N.  sin  a.  = k W.  + X.  sin  6.^,  - X.  sin  6. 
l ii  l c l l+l  i+1  l l 

+ Ei+1  cos  <5i+1  - E^  cos  6.^  (60) 

where  kc  is  the  critical  acceleration  required  to  bring  the  stresses 
on  the  probable  failure  surface  to  a state  of  limiting  equilibrium. 

The  limiting  equilibrium  condition  of  the  slip  surface  gives 

= (IT  - IL)  tan  <J>^  + c^k  sec  a..  (6l) 

where 

<t>!^  and  c!^  = shear  strength  parameters  at  the  base  of  the  slice  i 

U = force  exerted  by  the  pore  water  pressure  on  the  base 
of  the  slice 

b = width  of  the  slice 

Even  if  the  mass  contained  within  the  slip  surface  is  in  a state  of 
limiting  equilibrium,  the  mass  will  not  be  able  to  move  unless  shear 
surfaces  are  formed  within  the  sliding  mass.  Let  us  assume  that  the 
body  forces  X and  E on  the  slice  boundaries  are  such  that  they  are 
in  a state  of  limiting  equilibrium  as  well.  Then 

X ■■  (E  - PW)  tan  ?'  + c’d  (62) 
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where 

<t>'  and  c'  = the  average  strength  parameters  at  the  slice  edges 

FW  = the  force  exerted  by  the  pore  water  pressure  on  the 
side  of  the  slice 

Putting  Equations  6l  and  62  in  Equations  59  and  60  and  eliminating  IL  , 
* and  Xi+1  * we  obtain 


El*l  = *1  - Vo  * Vi 


where 


a.  = W.sinU’-c^+R.cos  ^+S.+1sin(^-ai-6i+1)-Sisin(^-a.-6i) 

i (gl*) 

cos  ( ♦ ! • - V* ! +1-«i+1 ) sec  ( * *+1 ) 


W.  cos  ( <t> ! - a.  ) 

X 1 1 


cos  (*•-«.+  ^+1  - 6.+1)  sec  *’+1 


cos  (<t>^  - sec  ^ 

cos  (♦•  - a±  + ?'+1  - «i+1)  sec  Vi+1 


R.^  = cp^  sec  ai  - tan  <J>! 


Si  = (cidi  ‘ Wi  tan  ♦P 


E=  a -bk  +Ee  = (a  +a  .e)-(b  +b  .e  )k 
n+1  n n c n n n n-1  n n n-1  n c 


+ E ,e  e i 
n-1  n n-1 

Therefore,  proceeding  further  with  the  recurrence  relation 


E = ( a +a  ,e  +a  „e  e ,,  + ...  to  n terms ) 
n+1  n n-1  n n-2  n n-1 


- k (b  + b e + b e e + ...  to  n terms) 
c n n+1  n n-2  n n+1 


+ E1  ' en  * Vl  • en-2  ei 


Since  in  the  absence  of  external  forces 


E,  = E =0 
1 n+1 


. a +a  ne  +a  0ee  ,+...to  n terms 
k = n n-1  n n-2  n n-1  

Q ~ ■ ■ ■ 

b +b  , e +b  _e  e , + . . . to  n terms 
n n-1  n n-2  n n-1 


Once  the  value  of  k is  determined,  all  E forces  can  be  determined 
from  the  recurrence  relation  (Equation  63)  and  then  Equation  62  gives  all 
the  X forces.  Equations  59  and  6l  will  give 

Ni  = (Wi  + Xi+1  COS  6i+l  ' Xi  COs  6i  - Ei+1  Sin  6i+l 

cos 

+ E sin  6.  - R.  sin  a.)  r (70) 

i ii  l cos  ($!  “ ai/ 

and  then  Equation  6l  will  give  all  values  of  T\  . 

53.  In  order  to  complete  the  solution,  we  have  to  obtain  the 
points  of  applications  of  N and  E forces.  For  this  we  make  the 
second  set  of  assumptions.  Let  us  assume  the  point  of  application  of 
all  E forces  is  known.  Then  taking  moment  about  the  corner  of  the 
slice,  we  obtain 

1i  = {Wi^Xgi  ~ xi^  " kcWi^ygi  " yi^  + Eizi  ' Ei+1  [zi+l  + bl  sec  “i  sin 
(“i  + 6i+l)]  + Xi+lbi  Sec  “i  COS  (“i  + ai+l)}/Nl  (71) 


where  (xg.^,  yg^)  are  the  coordinates  of  the  center  of  gravity  of  the 
slice. 

5U.  Alternatively,  we  may  determine  by  assuming  1^  as 

known,  in  which  case  the  last  1r  should  be  determined  as  part  of  the 
solution. 

55.  For  acceptability,  all  forces  (i.e.,  IT  , X.^  , and  ) must 

have  a positive  result,  and  1^  must  be  within  the  slice,  preferably 
within  the  middle  third. 

5 6.  In  obtaining  the  solution,  it  is  assumed  that  shearing  exists 
between  the  slices.  The  assumed  angle  6^  may  not  be  the  critical  one, 
and  it  is  therefore  necessary  to  try  various  angles  of  6^  . The  set  of 
6^  that  will  produce  minimum  kc  will  be  the  critical  set.  The  tech- 
nique of  achieving  this  minimum  value  will  be:  fix  all  but  one  6.^  , 
and  change  this  6^  gradually  to  obtain  the  minimum.  Then,  fixing 
this  angle  , proceed  to  another  point  until  all  the  6.^'s  are  tried 

and  the  smallest  value  of  k is  reached.  It  is  assumed  that  the 

c 

changing  of  the  order  of  6^  will  not  produce  a different  result.  How- 
ever, this  point  is  not  checked.  It  is  seen  that  the  assumption  of 
6^  = 0 for  all  slices  produces  a nonconservative  result,  but  only 
slightly.  Hence,  for  routine  procedure,  6.  = 0 can  be  used. 

57-  This  solution  can  be  converted  to  another  form:  fix  6^  = 0 
and  introduce  a local  factor  of  safety  FL  in  <f>'  and  c'  . It  is 
seen  that  as  F^  is  increased  from  1.0,  the  critical  acceleration  drops. 
The  minimum  k£  value  obtained  previously  corresponds  to  a local  factor 
of  safety  of  about  1.1  in  the  homogeneous  case.  However,  this  is  an 
observation  from  a very  small  number  of  test  cases  and  requires  further 
observation. 

58.  In  this  method  of  solution,  and  c'  , which  represent 

the  internal  strength  of  the  material  of  the  sliding  mass,  play  an 
important  part.  The  higher  the  values  of  and  c'  , the  higher  the 

critical  acceleration  is.  In  other  words,  the  higher  the  shear  strength 
of  the  material  inside  the  sliding  body  is,  the  lower  the  shear  strength 
required  at  the  base  of  the  slip  surface  to  maintain  static  limiting 
equilibrium.  Thus,  in  the  stability  analysis  of  the  existing  slip 


surfaces  where  the  strength  on  the  slip  surface  may  be  in  the  residuar 
state,  instead  of  in  the  peak  state  which  may  exist  inside  the  sliding 
body,  it  will  be  found  that  the  required  residual  strength  is  lower  than 
that  obtained  with  the  assumption  that  the  peak  strength  is  equal  to  the 
residual  strength.  The  effect  of  the  internal  strength  of  the  material 
could  not  be  studied  in  any  of  the  previous  stability  analysis  methods.* 

59.  Since  the  moment  equilibriums  of  the  slices  do  not  play  any 
part  in  the  determination  of  k£  (they  are  essential  for  the  complete 
solution),  the  slices  can  be  as  big  as  possible  and,  in  fact,  Eire  con- 
trolled by  the  change  of  the  angle  of  the  slip  surface.  Also,  the  posi- 
tion of  the  point  of  applications  of  E or  N forces  has  no  effect  on 
kc  , whereas  all  previous  methods  are  dependent  on  this  assumption.  The 
comparison  of  this  method  with  the  Sarma  (1973)  method  has  shown  that 
changing  the  point  of  applications  has  more  effect  on  the  production  of 
acceptable  results  than  on  the  value  of  kc  . 

60.  In  the  strictest  sense,  this  method  is  not  applicable  to  de- 
termining the  factor  of  safety  of  any  surface  except  the  probable  failure 
surface.  For  any  given  surface,  a factor  of  safety  exists,  but  the 
surface  may  not  fail.  However,  for  the  probable  failure  surface,  we  may 
use  the  factor  of  safety  to  mean  actual  reduction  of  the  strength  param- 
eters of  the  material  and  not  Just  the  strength  along  the  slip  surface. 

In  this  case,  the  method  is  applicable.  It  is  best  suited  for  existing 
slip  surfaces  where  kinematic  considerations  are  already  satisfied. 

61.  Figure  13  presents  an  example  of  this  method,  and  gives  its 
comparison  with  results  from  the  Sarma  (1973)  method  for  a homogeneous 
case. 

Method  2 

62.  This  method  is  given  here  to  show  the  ease  with  which  a rig- 
orous solution  of  the  stability  analysis  of  earth  dams  and  slopes  can  be 
obtained  in  terms  of  the  critical  acceleration  factor.  Therefore,  the 
method  will  not  be  developed  completely;  and  no  examples  will  be  given. 

* Sarma' s (1973)  method  determines  the  shape  of  the  X force  distribu- 
tion on  the  basis  of  the  internal  strength.  Because  of  the  factor 
X involved  in  the  solution,  the  exact  effect  is  not  deterministic. 
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Figure  13.  Comparison  of  method  1 with  Morgenstern-Price 

and  Sarma  methods 

63.  The  method  is  based  on  the  assumption  that  the  shape  of  the 
normal  force  distribution  on  the  slip  surface  is  known,  but  the  magni- 
tude is  not  known.  The  slip  surface  is  formed  of  segments  of  straight 
lines,  and  the  forces  on  the  base  of  the  slip  surface  are  as  shown  in 
Figure  lU.  From  the  vertical  and  horizontal  equilibrium  of  the  whole 
sliding  body. 


£(N!^  + IL)  cos  sin  (72) 

£T.  cos  a.  - £(NJ  + l^)  sin  = k,£w ± (73) 

The  Mohr-Coulomb  failure  criterion  gives 

Ti  = N1  tan  ^i  + cibi  sec  °i  (7*0 

U5 


Normal  stress 
distribution 


Figure  lU.  Stability  analysis  - method  2 

We  assume  that 

tq  = XPi  (T5) 

where  is  known  and  A is  to  be  determined.  Eliminating 

from  Equations  72,  73,  7^,  and  75»  we  obtain 

Y(W.  - c!b.  tan  a.  - U.  cos  a.) 

_ Lt  i i i _i l x 

cos  (<J>!^  - a.^)  sec  <J>!^ 

£(c!b.  - U.  sin  c^)  + A^P.^  sin  ( <|> !^  - c^)  sec  <J>^ 

“o 1 i -Ji 

6U.  Since  the  solution  depends  on  P.^  , the  following  expression 
for  P^  gives  acceptable  results  in  many  cases.  If  and  °^+1 

are  the  normal  stresses  on  the  two  edges  of  the  slip  line  segment,  then 

O!  + °<4.T 

Pj^  = - ^ sec  (78) 

where 

o'  = [cos  ♦*Yh(l  - Ru)  - c'  (sin  6 * sin  ♦ ')]  j-Tli^VsliMT  (79) 


N|  and 

(76) 

(77) 
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R 


u 


u 

yh 


(80) 


u = pore  pressure  at  the  point  under  consideration 

B = 2o  - (81) 

65.  In  order  to  compare  the  result  with  that  obtained  from  the 
other  method  of  slices,  it  is  necessary  to  divide  the  sliding  body  into 
vertical  slices  as  in  Figure  11  and  then  to  determine  the  interslice 
forces  and  their  points  of  application.  From  the  vertical  and  hori- 
zontal equilibrium  of  the  slice 


Xi+1  = (N£  + IL)  cos  cu  + TV  sin  or  - (82) 

E..,  = T.  cos  a.  - (N!  + U. ) sin  a.  - k W.  + E.  (83) 

l+l  1 1 i 1 1 ci  1 

Since  X is  known  from  the  solution,  IT  from  Equation  75  , and 

from  Equation  7**»  starting  from  known  initial  condition  of  = E^  = 0 , 

all  X and  E forces  can  be  calculated. 

66.  After  determining  the  X and  E forces,  the  points  of 
applications  of  E and  N'  forces  can  be  determined  by  assuming  the 
position  of  one  of  these  quantities  as  in  method  1. 

67.  This  method  appears  to  be  suitable  for  quick  computation  of 
k£  for  a slip  surface  and  should  give  acceptable  results  for  existing 
slip  surfaces.  However,  it  is  not  easy  to  see  whether  it  will  give 
acceptable  results  for  any  assumed  slip  surface.  Thus,  it  appears  that 
for  design  of  dams  where' one  has  to  start  with  assumed  failure  surfaces, 
Sarma's  method  (1973)  is  most  preferable  and  that  for  the  analysis  of 
existing  slips,  the  first  method  presented  herein  is  preferable. 


/ 


i 

< 

I 


Stability  of  Slopes  Under  Earthquake  Conditions 
with  Dynamic  Pore  Pressure 


68.  Under  the  earthquake  loading  conditions,  the  inertia  force 


4 


as  well  as  the  excess  pore  pressure  generated  by  the  earthquake  play  an 
important  part  in  the  stability  of  the  dam  and  its  foundations.  Until 
now,  this  excess  pore  pressure  was  introduced  into  the  analysis  as 
equivalent  static  pressure.  Seed  and  Martin  (1966)  proposed  using  the 
total  strength  determined  under  appropriate  cyclic  loading  conditions 
for  the  stability  analysis.  Seed  et  al.  (1969)  extended  this  idea  to 
develop  the  finite  element  method  for  studying  the  liquefaction  poten- 
tial of  the  dam  during  an  earthquake,  and  Seed  et  al.  (1973)  used  this 
method  for  the  analysis  of  the  Tower  San  Fernando  Dam  in  California, 
which  was  affected  by  the  San  Fernando  earthquake  of  1971.  However, 
since  this  is  a finite  element  method,  it  is  outside  the  scope  of  this 
report . 


69.  The  method  presented  herein  is  an  effective  stress  method  which 
can  take  into  account  the  excess  pore  pressures  developed  during  cyclic 
loading.  This  method  uses  the  concept  of  the  pore  pressure  parameters 
which  was  originally  defined  by  Skempton  (195*0  for  static  loading 
conditions. 

Dynamic  pore  pressure  parameter  Ar 

70.  Cyclic  loading  tests  on  saturated  samples  of  cohesionless 
soil  show  that  the  pore  pressures  increase  with  every  cycle  of  loading 
until  failure  by  large  strain  or  liquefaction  occurs  (Lee  and  Seed,  1967a; 
Lee  and  Seed,  1967b;  Peacock  and  Seed,  1968;  Martin  et  al. , 1975).  This 
pore  pressure  depends  primarily  on  the  applied  cyclic  shear  stress  and 
the  number  of  cycles.  Martin  et  al.  (1975)  proposed  a set  of  formulas 

to  determine  the  change  of  pore  pressure  ratio  Au/o^q  against  the 
number  of  cycles  for  a given  stress  ratio  At/o^o  in  a simple  shear 
test,  where  Au  is  the  pore  pressure,  is  the  effective  consoli- 

dation pressure,  and  At  is  the  applied  shear  stress.  The  same 
relationships  are  replotted  in  Figure  15  to  show  the  variation  of 
Au/o^o  against  At/o^q  for  a given  number  of  cycles.  From  these 
curves  we  see  that 


(8*0 
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Au/tTyo 


Figure  15.  Pore  pressure  response  in  simple  shear  test  (data 
from  Martin  et  al.,  1975) 

where  An  is  the  pore  pressure  parameter  that  is  a function  of  the  num- 
ber of  cycles  and  the  level  of  stress  ratio.  However,  for  the  initial 
part  of  the  curve,  the  curves  appear  linear;  and  An  is  a function  of 
the  number  of  cycles  only.  Higher  order  equations  for  these  curves  can 
easily  be  found  if  desired. 

71.  Initial  failure  is  defined  as  that  state  in  which  the  Mohr's 
circle  of  stress  first  touches  the  failure  envelope.  From  the  geometry 
of  the  Mohr's  circle.  Figure  16,  we  obtain 


where 


1 At_ 
R o' 


vo 


(1 


+ V 


(1  - 


sin  $'  - 


ko>‘ 


(85) 


(86) 


72.  This  relationship  is  shown  dotted  in  Figure  15  . Equations 


1+9 


J 


(assume  ^constant) 

Figure  16.  Initial  failure  in  cyclic  simple 
shear  test 

Equations  8U,  85,  and  86  are  derived  for  simple  shear  tests  where  the  ini- 
tial stress  conditions  are  represented  by  a'  and  k a'  in  which  k 

vo  o vo  o 

is  the  coefficient  of  earth  pressure  at  rest.  Equations  85  and  86  are 
derived  assuming  that  kQ  remains  unchanged  until  the  initial  failure. 

73.  In  the  slope  stability  problem  involving  dams  and  slopes, 
simple  shear  conditions  do  not  represent  the  true  state.  Neither  does 
the  triaxial  test  condition  where  either  no  rotation  of  the  principal 
stresses  are  allowed  or  they  are  rotated  through  90  deg.*  However, 
in  the  dam  during  earthquakes,  the  rotation  will  be  small  (Seed,  1968). 
Therefore,  the  consolidated,  undrained,  triaxiai  test  with  anisotropic 
consolidation  appears  to  be  nearer  the  truth,  and  we  may  assume  that 
the  pore  pressure  parameters  will  not  be  affected  by  the  rotation  of 
the  principal  stresses.  From  the  limited  test  data  available  in  the 
literature  on  triaxial  cyclic  test  results,  utilizing  both  isotropically 
and  anisotropically  consolidated  samples,  we  see  a similar  pattern  to 
the  simple  shear  tests;  and  we  may  write  for  the  initial  part 

* A table  of  factors  for  converting  U.  S.  customary  to  metric  (SI)  units 
of  measurement  is  presented  on  page  3. 
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and  for  the  anisotropic  consolidations,  perhaps 


f I 1 


3cfl+kc) 


’3c(l+kc) 


where 


o^c  = confining  cell  pressure  for  consolidation 
Ao^  = change  in  axial  stress 

* °ic^°3=  m) 

However,  Equation  88  will  need  extensive  laboratory  verification  to  see 


whether  or  not  An  is  independent  of  k£ 


The  parameters  in 


triaxial  tests  are  different  from  those  of  the  simple  shear  tests. 

7*t.  Equations  87  and  88  suggest  that  the  Skempton  ( 195*0  equa- 
tion can  be  adopted  for  dynamic  pore  pressure  as  well,  and  we  may  write 


Au  = B^Ao^  + An(Ao1  - Ao^)J 


where  B is  a pore  pressure  parameter. 

The  analysis 

75.  In  this  analysis,  the  suffix  o denotes  the  preearthquake 
conditions,  the  suffix  d for  dynamic  loading  conditions.  The  method 
is  based  on  Sanaa's  (1973)  technique  of  determining  kc  , the  critical 
acceleration  factor. 

76.  Let  us  consider  a possible  slip  surface  in  a slope  as  shown 

in  Figure  17-  We  subdivide  the  sliding  body  into  n vertical  slices. 

Before  the  earthquake  the  slice  is  acted  upon  by  forces  as  shown  in  the 

figure.  Let  us  assume  that  there  is  a factor  of  safety  Fq  under 

static  conditions  for  this  surface.  We  may  find  this  value  of  Fq  , 

the  interslice  body  forces  E and  X , and  the  normal  and  sheen: 

o o 

forces  at  the  base  of  the  slices  Nq  and  Tq  by  any  of  the  rigorous 
analysis  methods. 
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Aa  . 3>U/bsec6C 

* Mi/to  mc(C 
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Figure  17.  Stability  analysis  by  Sanaa  method 


77-  Let  us  now  assume  that  a horizontal  earthquake  load  kW  is 
applied,  under  the  influence  of  which  limiting  equilibrium  is  achieved. 
In  that  case,  k is  equal  to  kc  . Prom  Figure  17c,  we  may  then  write 


N.  cos  a + T ' sin  a = W - DX. 
da  d 


T.  cos  a - N.  sin  a = k W + DE. 
d d c d 


T = (N  - U - DU.)  tan  <J>'  + c'b  sec  a 
d d o d 


(91) 


(92) 


(93) 


where 

DU.  = the  dynamic  pore  pressure  developed  under  the  influence  of 

KW 

c 

W = weight  of  the  slice 
DX  and  DE  are  as  shown  in  Figure  17. 

78.  At  this  stage,  let  us  introduce  the  hypothesis  proposed  by 
Sarma  (1975).  This  hypothesis  states  that  if  a possible  failure  plane 
is  acted  upon  by  a normal  stress  o'  and  a shear  stress  x and  the 
shear  stress  has  a factor  of  safety  F , then  the  state  of  stress  in  all 
planes  will  be  the  same  as  that  in  which  the  material  has  a strength 
property  given  by  <|»'  = tan  "^(tan  4>'/F)  and  c'  = c'/F  • Under  this 
hypothesis,  it  is  possible  to  draw  the  Mohr's  circle  of  stress  for  static 
conditions.  Figure  17b.  From  the  geometry  of  the  circle. 


o'  = o'  + x (tan  <|>'  + sec  <»') 
X 9 0 o o o o 


3,0 


= o'  + x (tan  ' - sec  4>' 


(9«0 

(95) 


where 


c’  = c'/F 
° -1 

♦q  = tan  (tan  $'/F  ).  Adding  u (u  1 pore  pressure  at  the 
base  of  the  slice)  to  both  sides 
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(96) 


°l,0  ' 

°3,0  = 
a = 

Similarly,  for  the  dynamic 

°l,d  = 

°3,d  = 

We  define 

and 

and  therefore  Equation  90  gives 

Au  = b{a°3  + An(Ao1  - Aa^jJ 

Multiplying  both  sides  by  b sec  a , we  have,  with  ^ 

DUd  = B^(Nd  - Nq)  + Td|tan  4» ’ - sec  <J> ’ ( 1 - 2A)]  - TQ[tan 

- sec  ip ' ( 1 - 
o 

79.  Eliminating  Nd  and  Td  from  Equation  102,  with 
91  and  93,  we  obtain 

DU.  = a - b DX 
d d 

where 

- _ B(a  - Nq)  tdq- 
1 - BbQ  - eq 

51* 


a + t (tan  tp'  + sec  <p' ) 
o o o o 


o + t (tan  ip'  - sec  ip' ) 
o o o o 


stress  component 

loading  conditions,  we  shall  have 


od  + Td(tan  <Pd  + sec  ipd) 
od  + xd(tan  i Pd  - sec  <Pd) 


AO.  = ai  a ~ n 

1 l,d  1,0 


Ao  = O-  , - a 
3 3,d  3,0 


(97) 

(98) 

(99) 

(100) 

(101) 


Equations 

(103) 


(10U) 


Y 

I 


L 


f 


b = B(C  + fq)/l  - B(bQ  - eq)  (105) 

a - (W  + Uq  sin  a tan  ' - c'b  tan  a cos  <)>'  sec(<(>'  - a)  (106) 

b = tan  a tan  <J>'/(1  + tan  a tan  <J>')  (107) 

o 

C = cos  <J>'  sec(<(),  - a)  (108) 

d = (c'b  + W tan  $'  - Uq  cos  a tan  <J> * ) cos  $'  sec(<J>'  - a)  (109) 

e = tan  <f>'/(l  + tan  a tan  $')  (110) 

f = sin  4>'  sec(<|>'  - a)  (ill) 

q = tan  4>'  - sec  4>  * ( 1 - 2A^)  - (112) 

q = tan  <|»'  - sec  f (l  - 2A  ) (113) 

^o  o 0 n 


N = (W  + U sin  a tan  1 p'  - c'b  tan  a)  cos  U)'  sec  (^'  - a) 

00  00  00 

- DX  cos  ij>'  sec  (i|)'  - a)  (llU) 

00  o 

T = (c'b  + W tan  ii'  - U cos  a tan  ib ' ) cos  i|i'  sec  ( ' - a) 

0 0 00  ro  o To 


Now  eliminating  N,  and  T 
d 

DXd  tan  (<()'  - a)  + DE^  = 


- DX  sin  i|i'  sec  (♦'  - 0) 
o ro  o 


= tan 


1(t*n  6 


Co  F 


(115) 

(116) 

(117) 


from  Equations  91 » 92,  and  93,  we  obtain 


W tan  (<J>'  - a)  + (c'b  sec  a cos  «J> ' 


- Uq  sin  $')  sec  ($'  - a)  - (a  - bDXd)  sin  $'  sec  ($'  - a) 


and  rearranging  terms 

DX^  £tfan  ( 4> * - a)  - b sin  4> * sec  ( 4> * - a)]  + DEd  = D - kcW  (ll8) 

where 

D = W tan  (<t>’  - a)  + (c'b  sec  a cos  ’ - Uq  sin  <J>')  sec 

( 4> • - a)  - a sin  <t> ' sec  (<t>'  - a)  (119) 
summing  for  all  slices  and  since  £DEd  = 0 

£DXdQtan  (4> ’ - a)  - b sin  <f>'  sec  ( 4> ’ - a)}  + K^W  = £ D (12°) 

80.  Let  us  assume  that  the  point  of  application  of  Nd  forces 
are  known  and  that  the  coordinate  of  this  point  is  xm»ym  ^o:r 

practical  purposes , can  be  the  middle  point  of  the  base  of  the  slice. 
Then  taking  moment  about  the  centre  of  gravity  of  the  whole  sliding 
body  (x  ,y  ) and  remembering  that  the  net  moment  • W and  kW  about 
the  centre  of  gravity  is  zero,  we  have 

£(Nd  cos  a + Td  sin  a)(xm  - xg)  + £ (Td  cos  a 

- N.  sin  a)  (y  - y ) = 0 (121) 

a mg 

and  using  Equations  91.  92,  and  ll8  in  Equation  121, 

X“*d  { (lt.  - xg>  (yn  - V 

(122) 

The  problem  now  is  to  find  a set  of  Xd  forces  which  will  satisfy 
Equation  122. 

8i.  Let  us  assume  that  the  shape  of  the  Xd  force  distribution 


b sin  4> * sec  (<t>'  - a)]^ 


is  known  but  not  the  magnitude.  Then  we  can  say  that 


Xd=  XQ 


(123) 


where  Q is  known  but  X is  unknown.  Then  DX,  = X(Q_.  - Q. ) for  the 
th  ^ 1*^1  i 

i slice  = XP  Then  Equations  120  and  122  become 

X^P^tan  ( 4> * — ot ) — b sin  <j>'  sec  ($'  -a)]  + kc£w  =1 D (121*) 

x£p  {(xm  " xg)  + (yn  - yg)  Ctan  ' - “)  - b sin  sec  ( «(> * - a)[]j 

= Iw(xn  - xg)  + - yg)  (125) 


Therefore,  X and  k can  be  solved.  The  rest  of  the  solution  will 

c \ 

follow  the  same  pattern  as  in  Sarma  (1973)  and  is  not  shown  here. 
Procedure 

82.  A procedure  for  applying  the  above  method  of  analysis  can 
be  summarized  as  follows: 

a.  Determine  An  parameters  for  various  cycles,  say  1 
cycle,  5 cycles,  10  cycles,  20  cycles,  etc.,  from 
anisotropically  consolidated,  undrained  tests  for 
different  materials  of  the  dam. 

b.  For  each  set  of  An  values,  determine  the  k values 
for  the  possible  slip  surface.  Therefore,  we  obtain  the 
critical  accelerations  for  a corresponding  number  of 
cycles. 

£.  From  the  strong  motion  record  and  the  average  seismic 
coefficient  for  the  sliding  mass,  find  the  equivalent 
number  of  uniform  cyclic  seismic  coefficient  of  kc 
(Seed  et  al. , 1969).  If  the  available  number  of  cycles 
is  more  than  the  number  of  cycles  required  for  failure, 
the  factor  of  safety  will  drop  below  one.  Since  the 
available  number  of  cycles  will  depend  primarily  on  the 
choice  of  the  record,  it  is  better  to  consider  the 
critical  acceleration  for  10-20  cycles  as  the  most  crit- 
ical value  and  to  compute  the  probable  displacement  on 
the  basis  of  it. 
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PART  IV:  DISPLACEMENT  CRITERION  OF  DESIGN 


83.  Within  the  limit  equilibrium  principle,  a factor  of  safety 
smaller  than  one  represents  failure.  This  signifies  that  when  the 
factor  of  safety  is  less  than  one,  the  mean  strength  along  the  failure 
surface  is  less  than  that  required  to  maintain  equilibrium.  Therefore, 
a section  of  the  dam  or  embankment  will  slide  along  the  failure  surface 
and  will  come  to  rest  again  at  a time  and  place  when  the  new  mean 
stresses  do  not  exceed  the  strength.  Depending  on  the  amount  of  rela- 
tive displacement,  which  may  vary  from  a fraction  of  an  inch  to  a few 
yards,  the  structure  will  be  said  to  be  safe  or  to  have  failed.  It  is 
therefore  obvious  that  a factor  of  safety  less  than  one  cannot  be  per- 
mitted  under  the  static  conditions,  as  the  stresses  producing  this  stage 
will  exist  until  large  displacements  change  the  geometry  of  the  structure. 

81*.  However,  under  earthquake  conditions,  it  may  be  possible  to 
allow  the  factor  of  safety  to  drop  below  one,  as  this  state  will  exist 
only  for  a very  short  time.  Nevertheless,  the  consequences  of  allowing 
the  factor  of  safety  to  drop  below  one  must  be  known,  and  this  can  be 
measured  in  terms  of  relative  displacements  of  the  sliding  mass  relative 
to  the  main  body  of  the  structure.  Newmark  (1965)  proposed  the  sliding 
block  technique  to  measure  this  displacement,  later  Ambraseys  (1973) 
proposed  the  upper  bound  formula,  and  Sarma  (1975)  showed  that  displace- 
ments can  be  computed  on  the  basis  of  simple  pulses  once  the  critical 
acceleration  is  known. 

85.  Sarma' s (1975)  work  is  based  on  the  effective  stress  prin- 
ciple. The  displacements  were  computed  on  the  basis  of  the  assumption 
that  the  pore  pressures  at  the  base  of  the  sliding  block  remain  constant 
during  movement.  In  contrast  to  laboratory  testing,  where  pore  pres- 
sures continue  to  rise  after  the  factor  of  safety  has  become  equal  to 
one  in  cyclic  loading,  any  attempt  to  increase  the  earthquake  load  only 
increases  the  displacement,  creating  dilatency  in  the  slip  zone  and 
thereby  reducing  the  pore  pressure.  In  the  case  of  slope  failures, 
therefore,  the  assumption  that  pore  pressures  remain  constant  during 
movement  is  rather  conservative. 


* 

/ 
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86.  It  is  also  assumed  that  the  critical  slip  surfaces  can  he* 
converted  into  plane  slip  surfaces.  The  shear  strength  properties  of 
this  plane  surface  which  separates  the  sliding  block  from  the  main  body 
of  the  structure  are  the  average  and  c'  . The  pore  pressure  is 

equal  to  that  at  incipient  failure  (Figure  18). 


kg 

Figure  18.  Sliding  block  model 

87.  The  sliding  block  is  subjected  to  a horizontal  earthquake 
acceleration  kg  , which  imparts  a horizontal  load  kW  on  the  sliding 
block,  as  shown  in  the  figure.  The  earthquake  acceleration  is  greater 
than  the  critical  acceleration  kcg  , which  was  determined  for  the 
original  slip  surface. 

88.  The  driving  force  D on  the  sliding  block  is 

D = W sin  B + kW  cos  B (126) 

where 

W = the  weight  of  the  sliding  block 
B = the  slope  of  the  sliding  plane 
and  the  resisting  force  R is 

R = (W  cos  B - kW  sin  B - Uf)  tan  <t>'  + c'L  (127) 

where  L is  the  length  of  the  sliding  block.  When  the  earthquake 
acceleration  is  k£g  , the  driving  and  the  resisting  forces  are  equal, 
which  gives 
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W sin  B + k W cos  B = (W  cos  B - k W sin  B - U,)  tan  + c'L  (128) 
c cr 

where  Uf  is  the  force  due  to  pore  pressure  at  failure. 

Using  Equation  128  in  Equation  127 


R = W sin  B + k W cos  3 - (k  - k )W  sin  B tan  <t>'  (129) 

c c 


where  k is  the  seismic  coefficient  as  a function  of  time.  Since  the 
earthquake  acceleration  kg  is  greater  than  the  critical  acceleration, 
D is  therefore  greater  than  R . Under  the  action  of  the  net  driving 
force  (D-R),  the  block  will  move  down  the  slope;  and  by  Newton’s  Law, 
this  net  force  will  give  the  acceleration  X of  the  block  relative  to 
the  base.  Thus 


R 


= W ^ (k  - k ) (130) 

cos  <t>'  c' 

or 

X = 6 cos  ~ B)  (k  - k ) (131) 

cos  <t>'  c 

Since  the  earthquake  acceleration  kg  is  a function  of  time.  Equa- 
tion 131  can  be  evaluated  to  determine  the  total  displacement  x at 
the  end  of  the  earthquake. 

89.  For  a rectangular  pulse  of  duration  T/2  and  amplitude 

k g , Figure  19a,  the  solution  of  Equation  131  gives 
m • 


where 


Ux 


kmgT 


m cos  <fr'  1 / , ,,  . \ 

2 * cos  (<P'~  B)  2 km  kc 


(132) 


x = permanent  displacement  of  the  block  relative  to  the  base 
m 

k = maximum  value  of  the  seismic  coefficient 
m 

T = predominant  period  of  an  accelerogram  or  twice  the  duration 
of  the  pulse 
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Figure  19.  Forms  of  acceleration  pulses  used 
for  computation  of  displacement 


For  a triangular  pulse  of  duration  T/2  and  amplitude  k^g 
peak  at  time  XT/2  , Figure  19b, 


liy 

m cost’  _ 1 

k ^2  ' cos(<T  - 6)  = 2 S’ 
in 


k /k  )(1  - X k/k) 
cm  cm 


- x(k/k)‘ 

c m 


for 


0 - kc/km  - f1  ■ VcTttj] 


= |[(1  - k/kffl)3  (2+2V1  - X - 


for  (l-V^)/X  <kc/kffi<  1 

For  a half/sine  pulse  of  duration  T/2  and  amplitude  kmg 


, with  a 

f}  < W 

/x  (133: 

0 

(13k 

, Figure  19c 


12  ' cos'u^’e)  = (Kc/km  - Sin 


for  1 > k /k  > 0.725 
— cm 


= (k  /k  + a - it  + cos2a/2  cot  a/2)/irc 
V c m / 


(135) 


where 


0 1 kc/km  1 0.725 


q = a + kc/km  (cos  a - cos  q) 


(136) 


(137) 


a = sin-1(k  /k  ) 
c m 


(138) 


90.  The  solution  for  earthquake  records  is  obtained  by  assuming 
piece-wise  linear  acceleration.  Since  the  problem  is  that  of  a slope 
failure,  the  critical  acceleration  required  to  move  the  body  up  the 
slope  is  assumed  to  be  very  large,  and  no  reversal  of  displacements 
occurs. 

91.  Between  the  times  t..  , and  t.  , the  acceleration  k can 

i-1  1 g 

be  written  as 


k(ti  ■ ki-i  * 


(139) 


then  the  solution  of  Equation  131  gives 


T = (k,  ,-k  )(t-t.  ,)  + 


Vki-i  (t-Vir 


cosU'-fi)  g '“i-l  c"““i-l'  ’t^ 


(lUO) 


c°3,g  ■ . * = (k  _k  ) 

cosU'-B)  g ' i-l  c; 


(t-ti-i)  , Vki-1 


Wl 


+ ii_1(t-ti_1)  + xi_1 


(lUl) 


where  x is  the  dovnslope  velocity  of  the  sliding  block  relative  to  its 
original  position. 

92.  The  procedure  for  determining  the  displacement  for  earthquake 
record  will  then  be  to 

a.  Scan  the  earthquake  record  from  the  beginning  until  the 
pulse  is  bigger  than  k . This  will  give  the  initial 
conditions.  Determine  £he  velocity  and  displacement  at 
the  end  of  the  linear  pulse,  which  provides  the  initial 
conditions  for  the  next  pulse. 

b.  It  is  possible  that  the  movements  may  stop  in  the  middle 
of  a decreasing  pulse.  This  will  be  indicated  by  the 
change  of  the  sign  of  the  velocity;  in  which  case,  the 
time  for  zero  velocity  is  interpolated,  and  the  dis- 
placement is  computed  up  to  that  time  only. 

c_.  Once  the  movement  is  stopped,  it  will  not  start  unless 
the  pulse  is  bigger  than  k . Therefore,  we  are  back 
in  step  a. 


93-  The  total  displacement  at  the  end  of  the  record  gives  x 


m 


The  maximum  acceleration  in  the  record  gives  k 

____ — -m 


period  of  the  acceleration  spectrum  gives  T 

Ux  /k  gT2 
m m 


and  the  predominant 
We  can  therefore  compute 


The  earthquake  acceleration  time  history  is  not  symmetric 
about  the  zero  axis,  and  it  is  therefore  necessary  to  compute  the  dis- 
placement x by  changing  the  signs  of  the  accelerations. 

* o 

9h.  Figure  20  shows  the  quantity  (l/C)(l+x  /k  gT  ) computed  against 

m m 

k /k  for  the  four  cases  mentioned  above.  This  figure  shows  that  for 
cm 

k /k  > 0.5  , the  triangular  pulse  and  for  k /k  < 0.5  , the  rectangular 
cm  cm 

pulse  effectively  give  the  displacement  from  an  earthquake  record, 

where  k for  the  pulse  is  the  maximum  of  the  acceleration  record  and 
m 

T for  the  pulse  is  the  predominant  period  of  the  record.  Also,  it  is 

seen  that  the  absolute  displacement  for  k /k  >0.5  is  very  small  for 

c m 

real  earthquake  records.  C represents  cos  (<(>'  - 8) /cos  <J>  * . 

95-  The  use  of  the  sliding  block  technique  for  computing  the 
displacement  demands  the  following  from  the  dam: 

a. 

b.  An  average  o'  . 

c.  An  average  slope  angle  6 


A critical  acceleration  kcg 
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96.  The  methods  presented  in  Part  III  will  give  the  critical  ac- 
celeration kcg  . From  this  solution  we  can  draw  the  force-polygon  for 
the  entire  sliding  surface,  as  shown  in  Figure  21.  The  slope  o^  the 
shear  force  vector  can  be  used  as  the  slope  angle  3 . From  the  same 
solution,  we  can  obtain  the  average  shear  stress  t and  the  average 
effective  normal  stress  n'  on  the  slip  surface.  Since  we  can  compute 
c'  , the  average  angle  <J>'  can  be  obtained  from 


* 


av 


tan 


-1 


(1U2) 


5 = Shear  force 
R = Normal  n 


0 13)4  l.icuta" 


Figure  21.  Vector  diagram  of  forces 


PART  V:  ANALYSIS  OF  THE  LOPEZ  DAM 


97*  In  this  part,  the  Lopez  Dam  of  California  is  analyzed  by  the 
method  mentioned  in  the  previous  chapters,  and  the  stability  of  the  dam 
is  checked  against  the  San  Fernando  Earthquake  of  9 February  1971*  The 
probable  ground  motion  at  the  Lopez  damsite  during  this  earthquake  is 
supplied  by  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES). 
This  record  is  a modified  version  of  the  Pacoima  Dam  record.  The  record 
is  shown  in  Figure  22,  along  with  the  integrated  velocity  and  displacement 
records.  The  three  spectra  are  shown  in  Figure  23.  The  cross  section  of 
the  Lopez  Dam  is  shown  in  Figure  2k.  The  material  properties,  as  supplied 
by  the  WES,  are  shown  in  the  Table  2. 

Seismic  Coefficients 

98.  The  presence  of  the  debris  in  the  upstream  side  creates  an 
asymmetric  section.  However,  for  the  sake  of  analysis,  the  debris  will 
be  neglected  in  finding  the  inertia  forces  and  the  seismic  coefficients. 
The  average  shear  wave  velocity  for  a nonhomogeneous  section  can  be 
computed  in  the  following  way: 

b 

bl  b2  b3 

Sl’pl  S2’p2  S3’P3  ^ 

99-  Consider  an  elemental  slice  of  thickness  dy  . Assuming 
constant  strain  along  the  entire  width,  the  net  shear  force  T across 
the  section  is 


* <Gi»i  * 


°2b2  * 


°3V  W 


2 . 
iPibi 


(1U3) 


The  average  shear  stress  is  therefore 
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Figure  22.  Ground  motion  record  used  in  analysis  of  Lopez 
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Spectral  curves  for  ground  motion  record  of  Figure  22 


Figure  2h.  Lopez  Dam 


t 

t - — 
av  b 

and  therefore  the  average  shear  modulus  is 


(3 


G = ~ Ts2p.b. 

av  b **  1 1 1 

The  average  shear  wave  velocity  is  therefore 


However,  for  all  practical  purposes,  a linear  average  can  be  assumed 
so  that 


S 

av 


C 


Thus  for  the  Lopez  Dam  section,  considering  the  moist  densities,  the 


Table  2 

Lopez  Dam  - Material  Properties 


I 

I 


► 

r 


u 


Properties 

Zone  1 

Debris 

Zone  2 
Slope 

Zone  3 

Core 

Zone  4 
Foundation 

c*  , psf 

0 

0 

560 

224 

y 

4o.8° 

41.8° 

27.5° 

38.5° 

Ysatpcf 

130 

146.8 

139.3 

142.4 

Ymoistpcf 

107 

133 

123 

127.5 

S ft/sec 

700 

750 

850 

1000 

A * 

A1 

0.3 

0.3 

0.6 

0.2 

A * 

5 

1.0 

1.0 

2.0 

0.9 

A * 

10 

1.9 

1.9 

2.2 

1.7 

* 

O 

CVJ 

< 

2.2 

2.2 

2.5 

2.0 

* An  = pore 

pressxire  parameter 

A at  n 

cycles. 

average  shear  wave  velocity  is  795  ft/sec;  and  the  average  density  is 

o 

128  lb/ft  . The  shear  wave  velocity  of  the  foundation  layer  is  1000  ft/ 
sec,  and  the  density  is  127.5  lb/ft3.  Therefore,  m = S-jP^/S^  = 0.8 
and  q = m(H  - h^^/hp^  = 0.71  . For  this  pair  of  values  of  m and  q , 
the  first  root  of  the  transcendental  equation 


m tan  (qa  ) = — — — 

n J,  (a  ) 

1 n 

is  a = 1.34  , which  is  obtained  from  Figure  5 
puted  numerically  in  this  case). 

100.  Thus,  the  fundamental  period  of  the 


(148) 

of  Part  II  (also  com- 
dam -layer  system  is 


T = — = 0-27  sec  - (il*9) 

S a, 
av  1 

Figure  25  shows  the  seismic  coefficient  spectra  for  the  period  to  range 
from  0.2  to  0.4  sec  for  the  modified  version  of  the  Pacoima  Dam  record 
for  the  pair  of  values  of  m and  q mentioned  above. 

101.  If,  on  the  other  hand,  the  debris  is  not  neglected  but  is 


j 
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considered  as  part  of  the  foundation,  so  that  the  dam  height  is  reduced 
by  20  ft  and  the  layer  thickness  is  increased  by  the  same  amount,  then 
the  values  of  the  two  parameters  become  m = 1 and  q = 2.h  . For  this 
pair  of  values  a1  = O.5I+  , and  the  fundamental  period  of  the  dam-layer 
system  becomes 


= 0.37*  sec 


(150) 


102.  Figure  26  shows  the  seismic  coefficient  spectra  for  this  case 
for  the  period  to  range  from  0.2  to  G.U  sec.  From  Figures  25  and  26,  we 
can  replot  the  seismic  coefficients  for  the  two  periods  in  question  (Fig- 
ure 27).  The  maxima  of  the  two  will  give  the  design  seismic  coefficients. 
In  this  analysis,  this  is  obtained  entirely  from  the  first  case.  The 
point  at  the  crest  of  the  dam  is  obtained  from  the  point  response. 

r\  Crest  of  Dam 


Debris  level 


jC  40- 

O Layer  level 


Figure  27.  Seismic  coefficient 
versus  depth 


Lopez  Dam 

• case  I 

• case  2 
design 


05  I O 1-5  2 0 g 

Seismic  Coefficients 


* The  assumption  of  an  untruncated  wedge  is  not  strictly  applicable  in 
this  case.  However,  it  is  considered  to  be  a reasonable  approximation. 
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Critical  Accelerat ions 


103.  The  Lopez  Dam  section  is  analyzed  by  the  limit  equilibrium 
principle  using  the  Sarma  (1973)  method,  with  the  computer  program 
written  by  the  author  for  the  CDC  6L00-6600  computer  at  the  Imperial 
College  of  Science  and  Technology,  London.  The  various  probable  slip 
surfaces  and  their  critical  accelerations  are  shown  in  Figure  28. 

IOU.  WES  supplied  cyclic  triaxial  test  data  indicating  the  prob- 
able values  of  the  pore  pressure  parameter  for  various  cycles, 

shown  in  Figure  29;  and  the  results  are  tabulated  in  Table  2.  These 

A values  are  used  to  determine  the  critical  accelerations  for  a few 
n 

upstream  surfaces  wnich  were  considered  critical.  As  far  as  the  dy- 
namic pore  pressures  were  concerned,  only  the  material  below  the  water 
table  is  assumed  to  be  affected  by  the  earthquake  loading.  The  results 
of  the  analysis  are  also  tabulated  in  Table  3. 

105.  From  the  results  we  can  see  the  following: 

a.  There  will  not  be  any  postearthquake  slope  failure  due 
to  liquefaction.  Even  the  worst  case  can  withstand  a 
critical  acceleration  of  0.2  g with  20  cycles  of  pore 
pressure  development,  which  means,  in  other  words,  that 
without  an  earthquake,  there  will  be  a factor  of  safety 
greater  than  one  (a  critical  acceleration  of  0.2  g 
corresponds  to  a static  factor  of  safety  of  approxi- 
mately 1.6  without  earthquake). 

b.  Even  without  any  dynamic  pore  pressure  development,  the 
critical  accelerations  of  some  surfaces  are  smaller  than 
the  design  seismic  coefficients  (Figure  30).  Therefore, 
during  the  earthquake,  the  factor  of  safety  of  the 
structure  will  drop  below  one;  and  some  slip  surfaces 
will  develop. 


Displacements 


106.  Figure  31  shows  the  dimensionless  displacement  graph  for 
the  modified  version  of  the  Pacoima  Dam  record.  The  assumed  predominate 
period  is  0.U  sec.  There  are  two  peaks  in  the  acceleration  spectra, 
one  at  0.2  sec  and  the  other  at  0.U  sec;  however,  the  velocity  spectra 
shows  that  0.U  sec  is  perhaps  more  important  than  0.2  sec.  Table  U 


Cycles  of  stress 

b. 

Figure  29.  Probable  values  of  pore  pressure  parameter  A 


Table  3 

Results  of  Stability  Analysis 


Surface 

No. 

k 

c 

No  DPP* 

k 

c 

5 cycle 

k 

c 

10  cycle 

k 

c 

20  cycle 

6 

0.79 

0.59 

0.55 

0.53 

0.53 

7 

0.71 

0.51* 

O.Ul 

0.36 

0.35 

8 

0.66 

0.50 

0.3>+ 

0.28 

0.27 

9 

0.67 

0.U8 

0.33 

0.27 

0.26 

10 

0.85 

0.52 

0.37 

0.32 

0.30 

n 

0.70 

0.50 

0.U6 

0.39 

0.37 

12 

0.73 

0.51 

0.30 

0.2U 

0.22 

* Dynamic  pore  pressure. 
**  PP  = pore  pressure. 


Figure  30.  Comparison  of  criti- 
cal accelerations  and  design 
seismic  coefficients 


Table  k 


Computed  Displacements 


Displacement 


= ah 
ft 

k 

m 

Jv 

C 

No  DPP* 

K. 

C 

20  cycle 

c«* 

No  DPP* 
in. 

20  cycle 

e 

deg 

10 

1.25 

o.Uo 

— 

1.0 

7.2 

— 

20 

20 

1.15 

0.U5 

— 

1.0 

2.8 

— 

17 

30 

1.05 

0.66 

0.25 

1.0 

0.2 

Ih 

7 

Ho 

0.9U 

0.70 

0.30 

1.0 

0.06 

5.2 

2 

50 

0.80 

0.73 

0.22 

1.0 

— 

7.3 

0 

* Dynamic  pore  pressure. 

**  C = COS-^ — is  assumed  1.0. 
cos  4>  * 

shows  the  various  values  needed  for  computing  the  displacements  at  dif- 
ferent levels  and  also  gives  the  computed  displacements  from  Figure  31. 
From  the  table,  it  can  be  seen  that  the  nrjcimum  relative  displacement 
at  the  top  of  the  dam  will  be  about  7 in. , which  will  mean  a vertical 
drop  of  about  2.5  in.  and  horizontal  spread  of  about  6.5  in.  Deeper 
slides  will  have  negligible  displacement  without  the  dynamic  pore  pres- 
sure considerations. 

107-  With  a dynamic  pore  pressure  equi valent  to  20  cycles  of 
loading,  the  deeper  slip  surfaces  Eire  affected.  But  the  displacements 
are  still  smal.1,  the  maximum  being  about  lU  in.  This  displacement  does 
not  show  any  appreciable  drop  of  the  crest  level,  but  it  will  show  as  a 
crack  in  the  debris  at  about  100  ft  from  the  toe  of  the  dam. 

108.  Thus,  if  the  assumed  ground  motions  were  correct,  the  simpli- 
fied analysis  shows  that  the  Lopez  Dam  of  California  should  withstand  the 
San  Fernando  earthquake  of  9 February  1971  with  a very  minor  amount  of 
damage.  There  would  be  about  a 2.5-in.  slumping  of  the  crest  with  a cor- 
responding bulge  on  the  side  near  the  top  10  ft  of  the  dam.  There  should 
be  some  cracking  of  debris  about  100  ft  away  from  the  dam,  which  would  be 
caused  by  a small  movement  of  a deeper  slide. 


PART  VI:  CONCLUSIONS 


109.  The  method  of  analysis  and  design  of  earth  dams  under  earth- 
quake loading  conditions,  which  is  developed  in  this  report,  is  a simple 
one.  The  work  in  the  second  part  is  an  extension  of  the  original  work 
done  by  Ambraseys  and  Sarma  (1967),  which  now  includes  the  presence  of 

a foundation  layer.  The  presence  of  this  layer  does  seem  to  affect  the 
response  instead  of  merely  changing  the  fundamental  period.  The  rest  of 
the  conclusions  that  can  be  drawn  are  the  same  as  those  in  the  original 
paper  mentioned  above.  For  example,  the  upper  part  of  the  dam  is  af- 
fected more  than  the  lower  part,  and  a deeper  slide  has  a smaller  seismic 
coefficient  than  a slide  on  top  of  the  dam.  The  various  seismic  coef- 
ficient spectra  will  help  the  designer  in  computing  the  seismic  coef- 
ficients for  his  dam. 

110.  The  stability  analysis  method  with  the  dynamic  pore  pressure 
parameters  is  a simple  method  and  is  appropriate  for  design  purposes. 
However,  the  dynamic  pore  pressure  parameters  need  further  study.  Per- 
haps a high  order  curve  would  be  more  appropriate  than  a linear  one;  but 
with  the  present  available  data,  this  is  superfluous. 

111.  The  displacement  criterion  of  design  is  still  a matter  of 
opinion.  However,  the  displacements  computed  for  the  Lopez  Dam  do  not 
seem  to  disagree  with  the  field  behaviour. 

112.  The  analysis  of  the  Lopez  Dam  is  used  more  as  an  example  of 
the  concept  of  this  method  of  analysis  than  as  an  attempt  to  analyze  it 
in  detail.  For  example,  the  response  of  the  dam  and  the  seismic  coef- 
ficients are  not  shown  in  the  time  history  sequence;  only  the  maximum 
values  are  picked  up.  The  seismic  coefficients  for  particular  slip 
surfaces  are  not  computed  but  are  taken  from  appropriate  graphs.  Thus, 
the  analysis  shows  that  the  concept  is  a simple  one. 
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Applying  Laplace  transforms  to  Equations  A-l  and  A-2  with  the  conditions 
in  Equation  A-7, 


where  the  bar  signifies  the  Laplace  transform  of  the  function.  Let 


, and  SM-i 


(A-10) 


Equations  A-8  and  A-9  then  become 


, 1 d\  ,2—  — 

dy2  y dy"  Ul  6 


(A-ll) 


2—  2 

d »2  S12-- 
— 2 = - k u2  + g 

dy  s2 


(A-12) 


The  transformed  boundary  conditions  are 


dy 


= 0 at  y = 0 


(A-13) 


u2  = 0 at  y = H 


(A-lU) 


ux  = u2  at  y = h 


(A-15) 


du  du 

G —r—  = G — — at  y = h 
1 dy  2 dy  J 


(A-16) 


Solution  of  Equation  A-ll  with  Equation  A-13  gives 


ux  = CIQ(ky)  - ^2 
k 


(A-17) 


and  solution  of  Equation  A-12  with  Equation  A-lU  gives 


A3 


B sinh  — 


U2  ~ 


rsi  i /si  \ 

g-  k(y  - H)  cosh  I—  kyj 

(|“)  k2“’sh(|*H) 


(A-18) 


where  B and  C are  unknown  constants  and  I ( ) is  the  modified  Bessel 

o 

function  of  order  0 . 

Equation  A-15  then  gives 


sinh 


C = 


I^khT  • 


Pi  1 /si  \ 

<7=-  k(h  - H)  _ coshljpkHl 
■ ^ ^ Q ' ^ [_ 

m'j 


I (kh)  cosh 
o 


ft") 


Equation  A-l6  gives 


2 dUl  2 dU2 

SfPl  -r^  = S„p„  — t”  at  y = h 
11  dy  22  dy  J 


fl  !i 

S*  C2 


du^  du, 
dy 


_2 

dy 


Slpl 

c -rr  Ii(kh)  = B 
S2P2  1 


cosh 


[|t(—  H)]  ; 


where  I^( ) is  the  modified  Bessel  function  of  order  1. 


(A-19) 


(A-20) 


(A-21) 


Put 


Slpl 

S2p2 


= m 


(A-22) 


S-^h  - H) 


AU 


= a 


(A-23) 


(A-2U) 


I 


Equations  A-19  and  A-21  therefore  give 


— [" IQ(kh)  sinh  (bk)  - ml-Jkh)  cosh  (hk)"J 
kg  |_ml1(kh)  sinh  (ak)  - IQ(kh)  cosh  (ak)J 


Substituting  Equation  A-25  into  Equation  A-19  yields 


= Su-  . 


^2  * ml1(kh)  sinh  (ak)  - IQ(kh)  cosh  (ak) 
Putting  Equation  A-26  back  into  Equation  A-17  yields 


fc 


I (ky) 


2 I (h)  cosh  (ak)  - ml, (kh)  sinh  (ak) 
O X 


-] 


(A-25) 


(A-26) 


(A-27) 


Similarly,  substituting  Equation  A-25  into  Equation  A-l8  yields 

S S 

I (kh)  cosh  -r—  k(y  - h)  + ml.,  (kh)  sinh  -r—  k(y  - h) 
o Sg  1 bg 


- 1 


I (kh)  cosh  — k(h  - H)  - ml.,  (kh)  sinh  (ak) 
o b„  x 


(A-28) 


2.  The  inversion  of  Equation  A-27  and  Equation  A-28  can  be  ob- 
tained by  applying  Melin's  inversion  theorem.  We  can  write 


and 


U1  " « ’ F1 


u2  = g • Fg 


(A-29) 


(A-30) 


where  F1  and  F^  are  Laplace  transforms  of  F.^  and  Fg,  and 
therefore 


A5 


(A- 31) 


u ! = / \ - T>  dT 

° S1 


t 

u 2 = f 2 ^^2^  ~ T)  dT 


(A-32) 


Therefore,  the  problem  is  to  obtain  and  F^  from 


F =iL 

1 kHV 


kh)  cosh  (ak)  - ml1(kh)  sinh  (ak) 


F = — 
2 k2 


S S 

I (kh)  cosh  — k(y  - h)  + ml  (kh)  sinh  — k(y  - h) 

° S2  1 S2 

IQ(kh)  cosh  (ait)  - ml^kh)  sinh  (ak) 


(A-33) 


- 1 (A-3k) 


From  the  inversion  theorem 


S2  y+i»  ept  [ TQ(ry) 

F1  = 2^  (p2+cp)LIo<rh)  COSh  (ar)  ' mIl(rh)  Sinh  (ar] 


dp  (A-35) 


where  i = /-I  and 


S1  S1 

S2  Y+i*  pt  IQ(rh)cosh  g-  h(y-h)  + ml1(rh)sinh  g-  r(y-h) 
1 /*  6 2.  2 
2ni  (p2+cp)  Io(rh)cosh(ar)-mI1(rh)sinh(ar) 


where 


r = g-  (p2  + cp) 
S1 


- 1 dp  (A-36) 
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_____ 


To  apply  the  theory  of  residues,  the  poles  are  at  p=0;  p=-c;  and 


at  those  values  of  p = Pn  that  make 


I rh  cosh  ar  - ml,rh  sinh  ar  = 0 
o 1 


(A-37 ) 


The  roots  of  Equation  A-37  are  all  imaginary.  Let  rh  = ia^  he  the 


roots.  Then 


c + i/T 

Pn  ~ ~ 2 - 2 vc  " ^ ^ 


(A-38) 


are  the  poles,  and  Equation  A-37  becomes 


= m tan  qa 


(A- 39) 
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h Px  h 
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(A-Ul) 


c = 2X  u) 
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(A-k2) 


P = -X  id  + ico 
n n on  — on 


i1  - 


(A-U3) 


3.  The  residues  at  p = 0 and  p = -c  are  zero,  therefore 
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(A-UU) 
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APPENDIX  B: 
SEISMIC  COEFFICIENTS 


1.  For  the  one-parameter  sliding  wedge 

ah 

/y*n  dy 

- ^h — 

/ y dy 


rh  2Jo(any/h) 

J y a p (q,m,n)  ^ 
o no 


UJn(aa  ) 

1 n 

oa^P0(<l»m,n) 


(B-l) 


2.  For  the  four-parameter  sliding  wedge,  no  such  simple  form 
exists.  In  computing  <^Uq( ot ) , there  are  three  possibilities,  depending 
on  the  geometry  of  the  sliding  wedge. 

Case  1:  The  entire  sliding  surface  is  in  the  dam:  a 1 . 

Case  2:  Part  of  the  sliding  surface  is  in  the  foundation,  but 
the  exit  points  are  in  the  dam:  a > 1 . 

Case  3:  Part  of  the  sliding  surface  is  in  the  foundation,  and 
the  lower  exit  point  is  in  the  foundation:  a > 1 . 

3.  Sarma  (l 968)  determined  these  seismic  coefficients  by  comput- 
ing the  accelerations  in  the  dam  first  and  then  numerically  determining 
the  seismic  coefficient.  But  it  is  possible  to  formulate  the  seismic 
coefficients  analytically  for  the  three  cases  separately. 
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Case  3: 
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h ah 
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2 

Denominator  = ^B1  + 3^(1  - a^)  + B^(a  - l)J 


B1  a^h  B2  - B1  h 

Numerator  = — / yd>  dy  + — / y<|>  dy 
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APPENDIX  C: 

ROOTS  OF  THE  EQUATION 
m tan  (qan)  * J0(an^Ji^an^ 
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iu .... 


_ . _ '■ 
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ROOTS  OF  EQUATION  M.TAN(Q.AN)-JO(  AN)/Jl(  AN) 


M=  1.000 


0= 

.250 

.500 

.750 

1.000 

1.500 

2.000 

N 

1 

1.8774 

1.4934 

1.2238 

1.0312 

.7803 

.6262 

2 

4.3444 

3.6063 

3.1241 

2.7612 

2.2156 

1.8256 

3 

6.8703 

5.7708 

4.9329 

4.2832 

3.4266 

2.8824 

4 

9.4152 

7.8423 

6.7356 

5.8958 

4.7003 

3.8954 

5 

11.9457 

9.9259 

8.5349 

7.4405 

5.9713 

4.9731 

6 

14.4523 

12.0480 

10.3019 

9.0353 

7.2055 

6.0221 

7 

16.9515 

14.1309 

12.1211 

10.5883 

8.4874 

7.0509 

8 

19.4631 

16.2177 

13.9030 

12.1759 

9.T281 

8.1154 

9 

21.9868 

18.3294 

15-7036 

13.7332 

10.9907 

9.1633 

10 

24.5081 

20. 4l6l 

17.5025 

15.3172 

12.2530 

10.1980 

11 

27.0183 

22.5047 

19.2890 

16.8771 

13.4981 

11.2572 

12 

29.5235 

24.6117 

21.0967 

18.4584 

14.7684 

12.3407 

13 

32.0352 

26.7003 

22.8796 

28.0200 

16.0154 

13.3424 

14 

34.5547 

28.7901 

24.6856 

21.5997 

17.2758 

14.3988 

15 

37.0732 

30.8944 

26.4740 

23.1626 

18.5360 

15.4464 

16 

39.5845 

32.9841 

28.2716 

24.7411 

19.7849 

16.4858 

17 

42.0921 

35.0747 

30.0692 

26.3852 

21.0508 

17.5405 

18 

44.6041 

37.1771 

31.8591 

27.8827 

22.3005 

18. 5877 

19 

47.1218 

39.2677 

33.6619 

29.4473 

23.5597 

19.6286 

20 

49.6390 

41.3589 

35.4492 

31.0241 

24.8190 

20.6821 

M» 

.500 

Q- 

.125 

.250 

.375 

.500 

.750 

1.000 

n 

1 

2.2563 

2.1032 

1.9407 

1.7727 

1.4591 

1.2109 

2 

5.1492 

4.6135 

3.9859 

3.4996 

2.9465 

2.6526 

3 

7.9720 

6.7769 

6.0396 

5.6739 

5.1122 

4.4157 

4 

10.6397 

9.1594 

8.6100 

8.0771 

6.5901 

5.7817 

5 

13.1736 

11.8779 

11.0926 

9.8089 

8.5779 

7.5670 

6 

15.7793 

14.6492 

13.0064 

11.9487 

10.3526 

8.9194 

7 

18.5585 

17.2027 

15.2375 

14.3616 

12.0084 

10.7126 

8 

21.4531 

19.3642 

17.8545 

16.0986 

14.0852 

12.0594 

9 

24.3984 

21.7250 

20.1647 

18.2293 

15.5130 

13.8564 

10 

27.3519 

24.4396 

22.0479 

20.6448 

17.6754 

15.2000 

11 

30.2744 

27.2133 

24.4770 

22.3841 

19.1678 

16.9992 

12 

33.1114 

29.7715 

27.0719 

24.5112 

21.1383 

18.3406 

13 

35.7882 

31.9362 

29.1645 

26.9282 

22.9269 

20.1418 

14 

38.3206 

34.2913 

31.1715 

28.6689 

24.5716 

21.4821 

15 

40.9171 

37.0047 

33.7338 

30.7936 

26.6543 

23.2843 

16 

43.6907 

39.7788 

36.2293 

33.2120 

28.0797 

24.6234 

17 

46.5835 

42.3391 

38.1483 

34.9530 

30.2411 

26.4261 

18 

49.5287 

44.5038 

40.3674 

37.0763 

31.7365 

27.7644 

19 

52.4821 

46.8576 

42.9844 

39-4948 

33.7031 

29.5684 

20 

55.4064 

49.5701 

45.3016 

41.2369 

35.4950 

30.9058 

M=  0 


0-  0 

1 2.4048 

2 5.5201 

3 8.6538 

4 11.T915 

5 14.9309 

6 18.0710 

7 21.2163 

8 24.3525 

9 27.4933 

10  30.6346 

11  33.7762 

12  36.9171 

13  40.0584 

14  43.1998 

15  46.3412 

16  49.4820 

17  52.6241 

18  55.7657 

19  58.9070 

20  62.0485  C2 


APPENDIX  D: 

ACCELERATION  RECORDS  OF  FIVE  EARTHQUAKES 


* 


In  each  of  the  following  figures  (D1-D9),  a represents  accelera- 
tion, velocity,  displacement  records,  and  energy  flux  plots;  and  b rep- 
resents acceleration,  velocity,  and  displacement  spectra  of  strong  motion 
records . 
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APPENDIX  E: 
SEISMIC  COEFFICIENTS  AND  POINT 


ACCELERATION  SPECTRA 
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APPENDIX  G: 
NOTATION 


PART  II 


n 

b 

B 
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max 

F1 
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g(t) 
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The  n 
th 


th 


root  of  the  equation  J (a  ) = 0 

on 


o 
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n ~ root  of  the  equation  m tan  q(a  )/J^(a  ) 

Width  of  crest 
Unknown  constant 

Damping  coefficient  for  both  dam  and  layer 

Unknown  constant 

Total  inertia  force 

Maximum  inertia  force 

A function  defined  in  the  text 

A function  defined  in  the  text 

Acceleration  due  to  gravity 

Base  accelerations 

Maximum  base  acceleration 

Shear  modulus  of  the  material  of  the  dam 

Shear  modulus  of  the  material  of  the  layer 

Height  of  dam 

Total  height  from  crest  of  dam  to  rock  level 
Square  root  of  -1 

Modified  Bessel  function  of  order  zero 

Modified  Bessel  function  of  order  one 

Bessel  function  of  the  first  kind  of  order  zero 

Bessel  function  of  the  first  kind  of  order  one 

Average  seismic  coefficient  expressed  as  a fraction  of  gravity 

Average  seismic  coefficient  expressed  as  a fraction  of  the 
maximum  ground  acceleration 

Length  of  dam 

The  ratio  Slpl/S2p2 

Total  mass  of  the  sliding  body 

A function  defined  in  the  text 

Amplification  factor,  magnitude  of  ground  acceleration 
A function  defined  in  the  text 
The  ratio  mfH-lOpg/hp^) 


G2 


1 


w 

y 

a 

X 


+n(y) 

*V>\, 


*n<y> 

<u  ,w 
n on 


Response  acceleration 

Shear  wave  velocity  in  the  dam  and  in  the  layer,  respectively 
Time 

Horizontal  displacement  of  a point  in  the  dam  and  in  the 
layer,  respectively,  relative  to  the  rock  base 

Velocity  of  a point  in  the  dam  and  in  the  layer,  respectively, 
relative  to  the  rock  base 

Absolute  acceleration  of  a point  in  the  dam  and  in  the  layer, 
respectively 

Weight  of  the  sliding  mass 
Vertical  coordinate 

The  depth  of  the  base  of  the  sliding  surface  from  the  crest 
Damping  as  a factor  of  critical 
Mass  density 

Mass  density  in  the  dam  and  in  the  layer,  respectively 
Angle  of  slope  with  the  horizontal 
n mode  shape  for  the  dam 

Shape  function  for  one-parameter  and  four-parameter  sliding 
wedge,  respectively 

n mode  shape  for  the  layer 

Damped  and  undamped  circular  frequency,  respectively,  of  the 
n^*1  mode  of  the  dam- foundation  system 


PART  III 

A Pore  pressure  parameter 
An  Dynamic  pore  pressure  parameter 

b Width  of  the  slice 

B Pore  pressure  parameter 

c'  Cohesion  at  the  base  of  the  slice  in  terms  of  effective  stresses 
c'  Average  cohesion  at  the  side  of  the  slice 
c|  Cohesion  at  the  base  of  the  slice  i 
DU^  Dynamic  pore  pressure  developed  under  the  influence  of  kcW 

E Total  normal  force  on  the  side  of  the  slice 

Eq  Inter slice  body  force 

F Factor  of  safety 


Part  III  (Continued) 


x »y 
g g 

x ,y 
m Jm 


Static  factor  of  safety 

Local  factor  of  safety  in  the  side  of  the  slice 

Height  of  the  section  above  the  slip  surface 

Slice  number 

Seismic  coefficient 

Horizontal  earthquake  load 

Critical  acceleration 

Coefficient  of  earth  pressure  at  rest 

Coefficient  of  anisotropic  consolidation 

Distance  of  the  point  of  application  of  the  normal  force  at  the 
base  of  the  slice  from  the  edge  of  the  slice 

Normal  force  at  the  base  of  the  slice  in  terms  of  total 
stresses  and  effective  stresses,  respectively 

Normal  force  at  the  base  of  the  slice 

Force  exerted  by  the  pore  water  pressure  on  the  side  of  the 
slice 

Shear  force  at  the  base  of  the  slice 
Shear  force  at  the  base  of  the  slice 
Pore  pressure  at  the  base  of  the  slice 

Force  exerted  by  the  pore  water  pressure  on  the  base  of  the 
slice 

Weight  of  the  slice 
Horizontal,  coordinate 

Coordinate  of  the  center  of  gravity  of  the  whole  mass 
Coordinate  of  the  point  of  application  of  the  normal  force 
Shear  force  on  the  side  of  the  slice  section 
Interslice  body  force 

Height  of  the  point  of  application  of  E force  above  the 
slip  surface 

Slope  of  the  base  of  the  slice  with  respect  to  the  horizontal 
Density  of  the  material 

Slope  of  the  base  of  the  slice  with  lespect  to  the  horizontal 
Change  of  pore  pressure  due  to  cyclic  loading 


Part  III  (Continued) 


Ao^  Change  in  axial  stress 

Ao^  Change  in  cell  pressure 

At  Applied  cyclic  load  in  simple  shear  test 
X An  unknown  to  be  determined  from  the  solution 

a Stress  component 

o'  Normal  stress  at  the  base  of  the  slice 
Vertical  consolidation  pressure 
o|c  Axial  pressure  for  consolidaton 

o^c  Confining  cell  pressure  for  consolidation 

t Shear  stress 

<f>!  Shear  strength  parameter  at  the  base  of  the  slice  i 
<t>'  Average  strength  parameter  at  the  edge  of  the  slice 
PART  IV 

c'  Average  cohesion  in  terms  of  effective  stresses 

av 

D Driving  force 

k The  maximum  value  of  the  seismic  coefficient 

m 

L Length  of  sliding  block 

n'  Average  effective  normal  stress  on  the  slip  surface 
R Resisting  force 

T Predominant  period  of  an  accelerogram  or  twice  the  duration 
of  the  pulse 

Force  due  to  pore  pressure  at  failure 
W Weight  of  the  sliding  block 

x Instantaneous  displacement  of  the  sliding  block  relative  to 
the  base 

x Instantaneous  velocity  of  the  sliding  block  relative  to  the 
base 

X Instantaneous  acceleration  of  the  sliding  block  relative  to 
the  base 

x Permanent  displacement  of  the  block  relative  to  the  base 
m 

0 Slope  of  the  sliding  plane 

T Average  shear  stress  on  the  slip  surface 
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